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1. INTRODUCTION 


THE branch of mathematical physics known as the theory of the elasticity 
of solids is based on certain notions regarding elastic stresses and strains 
which have been accepted doctrine for many years. It is therefore rather 
surprising to find that those notions are not sustainable and that the theory 
based on them has to be laid aside in favour of a stricter formulation. The 
need for such revision was broadly indicated in a recent publication in these 
Proceedings by the present authors (1955); but in view of the importance 
of the matter, the present paper is prefaced by a clearer and more precise 
exposition of the necessity for a revision of the classical theory. ° 


Our recent paper dealt with the case of isotropic solids in a formal 
manner, considering them as homogeneous substances exhibiting the highest 
possible symmetry in their elastic properties. Usually, however, the so- 
called isotropic solids are merely polycrystalline aggregates, and a discus- 
sion of their elastic behaviour should therefore properly be based on a 
_. consideration of the nature and properties of such aggregates. This is a 
task which we hope to be able to address ourselves in the not-too-distant 
future. In the present paper we shall consider the case of truly homogene- 
ous but anisotropic solids, in other words, crystals. The subject will be 
dealt with from the phenomenological standpoint since this proves to be 
entirely adequate. We may remark, however, that the conclusions reached 
are in complete accord with the results of the atomistic approach to the 
theory of elasticity as developed both from the static and dynamic stand- 
points in a paper by one of us (Viswanathan, 1954) and illustrated by a 
detailed discussion of a particular case of great interest, viz., diamond. 


2. THe PHysicAL CHARACTERS OF ELASTIC STRESS AND STRAIN 


A solid body may be defined as a material in which the volume elements 
fetain their relative positions and orientations alike when the solid is at 
test and when it is in a state of movement, translatory or rotatory, as the 
case may be. Elastic stresses and strains arise when the situation thus de- 
scribed is departed from to ever so small an extent, Hence, it is evident 
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that the theory of elasticity is concerned with changes in the relative positions 
and/or the relative orientations of the volume elements of the solid, in other 
words, with linear displacements and/or. angular movements of the volume 
elements with respect to their neighbours. This is illustrated by the familiar 
examples of a straight bar which is stretched or bent or twisted by the applica- 
tion of appropriate external stresses. From the general considerations 
indicated above as well as from the particular examples mentioned, it is 
clear that only in very special cases would elastic strains be such that a 
straight line drawn through the solid in any direction in the unstrained condi- 
tion remains straight in the strained state, suffering only elongations or con- 
tractions. In the general case, and inevitably so when the strains involve 
differential rotations, a straight line in the unstrained condition would be 
curved in the strained state, and such curvature cannot possibly be ignored 
in the theory. Hence, it is clear that we have, in general, to consider strains 
and stresses which are heterogeneous, in other words, strains and stresses 
whose specifications vary from point to point within the solid. These varia- 
tions necessarily enter into the equations of equilibrium in the static state 
and into the equations of motion in dynamic. behaviour. 


3. ANALYTICAL SPECIFICATION OF STRESSES AND STRAINS 


The mathematical theory of elasticity proceeds on the basis that the 
strains and stresses in the interior of the solid can be expressed in terms of 
the movements of the smallest possible elements of volume into which it can 
be imagined to be subdivided and of the forces acting on them. If the volume 
elements be small enough, their movements can be described completely 
in terms of the three positional co-ordinates of each element and their varia- 
tions. Likewise, when the elements of volume are small enough, the inter- 
actions between each element and its neighbours can be expressed in terms 
of tractive forces alone, it being then clearly unnecessary to introduce any- 
thing in the nature of couples or torques. On the basis of these ideas, the 
state of strain in the solid at any given point can be expressed by resolving 
the displacement of the elementary volume originally located at such point 
along three mutually perpendicular directions and differentiating these three 
components of displacement again along each of the three axes in turn. We 
thus obtain the nine components of the strain tensor. Likewise, for speci- 
fying the state of stress to which the volume element is subject, we consider 
the tractive forces acting on an infinitesimal area drawn respectively normal 
to the three co-ordinate planes in turn at the position of the element and 
then again resolve these tractive forces along each of the three co-ordinate 
axes, We thus obtain the nine components of the stress tensor, 
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It is evident that this method of representation of the stresses and strains 
uniquely defines the state of the solid at any point and also enables us to 
determine whether the element of volume would or would not remain in 
equilibrium. If the components of stress do not vary along any of the three 
axes, the element would necessarily remain at rest. If, on the other hand, 
the stress components vary, their differential coefficients along the normals 
to the planes on which they act gives us a measure of the forces on the volume 
element in their respective directions. Adding up the three forces along 
each axis thus evaluated and putting their sums separately equal to zero, 
we obtain the conditions of equilibrium. 


The well-known and familiar treatments given in the standard treatises 
proceed on the basis that the components of the strain and stress tensors 
are both reducible in number from nine to six. The arguments on which 
the reduction in number from nine to six of the components of strain is 
based may be summarised by the statement that the elastic strains can be 
separated into what are called “ pure strains” and “ rotations,” and that 
the latter can be ignored. That this argument is unsustainable will be evi- 
dent at once from the remarks made in the foregoing section regarding the 
physical nature of elastic strains. We have, in general, to take account of 
both differential displacements and differential rotations and it is therefore 
not permissible to eliminate the rotational parts of the strain, these being 
physically quite as real as the irrotational parts. 


The arguments justifying the reduction in the number of the inde- 
pendent components of stress from nine to six are based upon the idea that 
equilibrium would be possible only if the angular momenta of the tractions 
taken about each of the three co-ordinate axes in turn cancel each other out. 
That this idea is misconceived will be evident from the remarks already 
made earlier regarding the conditions necessary for equilibrium. In the 
case of homogeneous strains, the tractive forces acting on each volume ele- 
ment necessarily balance each other. In the case of heterogeneous strains, 
the conditions of equilibrium can be expressed in terms of the differential 
coefficients of the stress components along the normals to the planes on which 
they act, as already explained. In either case, if the equilibrium conditions 
for each volume element of the solid are satisfied, then the solid as a whole 
necessarily remains in equilibrium; vice versa, if the external stresses acting 
on the solid are such that it remains as a whole in equilibrium, the elastic 
stresses would everywhere necessarily be such as to ensure equilibrium of 
the individual volume elements. It follows that no general relations con- 
necting the magnitude of the tensor components and enabling their number 
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to be reduced from nine to six can be derived from considerations based on 
the conditions for equilibrium. 


We may summarise our conclusions by stating that neither the reduc- 
tion of the strain components nor the reduction of the stress components 
in number from nine to six has any theoretical justification; a correct and 
complete theory of elasticity has necessarily to take all the nine components 
of the stress and strain tensors into consideration. 


4. THe STRESS-STRAIN RELATIONSHIPS 
Writing the nine components of the strain tensor as 
Ure Uyy Uzz Uyz Uzy Uzx Uxz Uxy Uyx 
U, Us Us Ug Us Ug Uy Ug Uy 
and likewise the nine components of the stress tensor as 
Tox Tyy Tzz Tyz Tzy Tzz Txz Txy Tyx 
7 Se Se te le he Ba Be 


the stress-strain relations can be expressed in the general form 


=_—— — (1) 


n=1 


and involve 8/ constants. Here the constant dm,» relates the stress Tm to 


the strain wu, and is the ratio of the two for a deformation in which all 
components other than uy vanish. 


The 8/ constants figuring in (1) are not all independent, but reduce in 
the first instance to forty-five for all solids in view of the relations 


dmn = Inm (m, n= 1, 2....9) (2) 


The above relations follow from the well-known theorem of reciprocity 
relating forces and the corresponding displacements of dynamical systems. 
The reciprocity relations further enable us to write down the expression for 


the deformation energy per unit volume in the neighbourhood of any point 
and this is given by 
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Thus in the general case of a completely anisotropic solid, we have 
forty-five elastic constants instead of the 2] contemplated by the classical 
theory. 

TABLE I 





No. of Elastic 


Constants 
Crystal System Point Group 





I II 





Triclinic All 21 45 





Monoclinic All 25 





Orthorhombic All 15 


C, (4) 
S, (4) 
Caa(4/m) 





Tetragonal 


Ca.(4mm) 
Sa=Deoe (42m) 
Dy, (4/m2/m2/m) 


C; (3) \ 
Se (3) 


D, (32) 


D, (422) 








Trigonal 


Cy, (Grn) 
Das (32/m) 


C, (6) 
Cy, (6) 
Ce, (6/m) 


D, (622) 
Coe (6mm) 











Hexagonal 
Dg, (62) 
Dg, (6/m2/m2/m) 


T (23) _ 
T, (2/m3) 


O (432) 





Ts (43m) 
O, (4/m32/m) 
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The number of independent elastic constants which is forty-five for a 
completely anisotropic solid diminishes in the case of solids possessing the 
various elements of symmetry characteristic of the different crystal classes, 
coming down to four for crystals of the classes Tg, O and Op. The number 
of surviving constants in each symmetry class.can be computed in an ele- 
gant fashion by adopting the group-theoretical method developed by 
Bhagavantam (1949). One is concerned in the present case with a linear 
relationship between nine stress components and nine strain components, 
the constants of proportionality being the elastic constants (matrix dm,). 
Further, the elements of the “ elastic constants 9 by 9 matrix” satisfy the 
relation dmn=dnm-. With these restrictions, the transformation matrix 
for the elastic constants can be written out and the corresponding character 
can be deduced. This comes out as 


X;’ (R) = 1léc* + 16c® + 8c? + 4c + 1 (4) 


where c = cos ¢, R is asymmetry operation and ¢, the rotation. The plus 
sign is used for proper rotations, and minus for an improper one. The 
corresponding character for the 2l-constant theory is 


Xj’ (R) = 16c4 + 8c — 4c? + 1 (>) 


Table I gives the number of independent constants according to the two 
formule for each symmetry class, while the constants that survive and those 
that vanish are exhibited in detail for the various cases in Tables II to XII. 
A comparison of these tables amongst themselves will enable the reader 
to realise how the existence of common symmetry elements results in the 
appearance of common features in the Tables of elastic constants. For 
example, all the twenty constants that vanish for monoclinic crystals also 


TABLE II 





| Triclinic (C,, C;) | 








| dy dys dis dy, dys dy, dy, dis diy 
| 22 dos doq dos 26 do, dog dog 
33 ds4 35 ds dz sg sg 
dy, dis dag day das ag 
55 56 ds. dsg dso 

66 dey deg deg 

dy, dag drg 

45 constants deg deg 











J 
<< 
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TABLE III 





Monoclinic (C,, C2, Cs,)—diad axis || z, plane of reflection + to 





















































dy, ays dys 0 0 0 0 dg dg 
do» dog 0 0 0 0 dog dog 

a 0 0 0 0 dsg deg 

Aaa das Ay6 dy, 0 0 

dss 56 ds, 0 0 

66 der 0 0 

daz 0 0 

25 constants dsg deg 
dog 

TABLE IV 
Orthorhombic (C.,, De, Dox) 

dy dys dis 0 0 0 0 0 0 

es dog 0 0 0 0 0 0 

ds 0 0 0 0 0 0 

dys das 0 0 0 0 

dss 0 0 0 0 

dee de, 0 0 

daz 0 0 

15 constants dsg deg 
yg 

TABLE V 
Tetragonal (C,, S4, C4,)—tetrad axis || z 

dy ds, dys 0 0 0 0 As dy 
dy 13 0 0 0 0 —dyy —dys 
dog 0 0 0 0 ds —Adsg 

daa das dag 0 0 0 

dss 0 —Adx 0 0 

ds. das 0 0 

daq 0 0 

ds 











13 constants 
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TABLE VI 





Tetragonal (Cy, Sq, Dy, D4,)—tetrad axis || z 











0 
0 
0 
0 
0 
Aas 





9 constants 








f eeeeooe 


nRNoco°cece 


~ 
& 





TABLE VII 





Trigonal (Cs, S,)—triad axis || z 








dg = —d12— Ago das 


15 constants 


dys Aye dy, 
dy dye = —dy 
0 
daa dy, 0 
45 
sg 








TABLE VIII 





Trigonal (Ds, Dsa, Csp)—triad axis || z 








dy dy, 
d, 1 


deg=4y —d.—Agg 


dys dy, 


dis ed | 


0 
0 
dss 0 
dus : 

das 





10 constants 


55 
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TABLE IX 





Hexagonal (C3, Cg, Ce,)—hexad axis || z 





‘SS 
eo 


0 0 
0 0 
0 0 
das d 
dss 


— 
ao 


eS 


46 
0 —dyg 
das ds 


55 








sooCoo 


ao 


11 constants 














TABLE X 





Hexagonal (Dg, Cer, Dg, Dg,)—hexad axis || z 





dy dy, dis 
dy, 13 
33 





deg=dyy —dy2—Agg 


0 
0 
0 
0 
0 
dys 
dss 


8 constants 


goooseeso 
AxRcooocooce 














TABLE XI 





Cubic (T, Ty) 








Nococooococeo 


— 
ss 


5 constants 





eBoooooeseo 


a 
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TABLE XII 





Cubic (Tz, O, O,) 


d, 1 d, 2 d, 2 








4) 0 0 0 0 0 
u dye 0 0 0 0 0 0 
dy 0 0 0 0 0 0 
dig das 0 0 0 0 
ds 0 0 0 0 
dys das 0 0 
dag 0 0 
4 constants das das 
d 











— 
- 











disappear in the orthorhombic, tetragonal and cubic systems. Likewise, 
all the thirty constants that vanish for orthorhombic crystals vanish also 
for those tetragonal crystals which possess three mutually perpendicular 
diad axes and for all cubic crystals. A noteworthy feature is that the cubic 
crystals which exhibit a four-fold axis of symmetry have only four different 
elastic constants, whereas those cubic crystals that do not exhibit this feature 
have five constants different from each other. 


5. 'WAVE-PROPAGATION IN CRYSTALS 


The general equations of motion of an elastic solid are given by 


aS * iy 
uy Tyr , Tyy , Wyz 
‘Fr “a = Se (6) 
en * ee 
where p is the density of the medium and T;,,.... are the stress components. 


If the solid is in equilibrium, the quantities on the left-hand side of these 
equations, representing the acceleration of an unit volume element, will vanish 
and we get the conditions of equilibrium of the solid. 


To evaluate the velocity of propagation of waves in the solid in any 
given direction, we first seek solutions of the above equations which are in 
the form of plane waves of the type 


u = Aexpi = (vt — e.r) (7) 
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Denoting the components of A in the directions of the co-ordinate axes by 
A*, AY, A? and using once again the four-suffixed symbols for the elastic 
constants, we now get on substitution of (7) in (6) that 


pv? A* = XS dz, yy ez ey AY (8) 
x yy 


and two similar equations for the y and z components. 


For a wave progressing in the direction (/mn) equations (8) can alter- 
natively be written as 


(Ara — pv?) A® + Ay AY + Azz A? = 0 


Any AZ ++ (Ayy — pv?) AY + Ayz AZ = 0 (9) 
Azz A® + Ayz AY + Azz A? = 0 
where Axx, Ary.... are given by the scheme 
Axx dy dsg dry 2dr 2d1, 2d; P 
Ayy Agy don gy 2do4 2dag 2dog m* 
Azz sa deg 5 dss 2ds 2dsg 2ds6 n* 
Ayz des dos dy4 (dog + das) (dg + di) (dsg + dog) mn 
Azx Aig dss sz (dsz + dg) (der + dys) (dis + deg) nl 
Ary sy deg day (doz + dys) (dig + dry) (di, + dg9) lm 


Equations (9) determine the velocities of propagation of the three types 
of waves in any direction for crystals of the triclinic system which possess 
no symmetry of structure at all. The number of constants figuring in the 
wave equations will diminish rapidly as we pass on to crystals of higher 
symmetry, and become only three for crystals of the Tg, O and Op» classes. 
We tabulate below the wave equations for the different classes of crystals 
taking into account of their symmetry. 


I. Monoclinic system (Cs, C2, Cp)—diad axis parallel to the z-axis 
(dil? + deg m® + dzqn® + 2dg lm — pv’) A® + {dygl? +- dagm® + dyyn* 
+ (diz + deg) lm} AY + {(dsg + a5) mn + (ds + de) nl} A? = 0. 
{dyol? + dag? + dygn® + (diz + dag) lm} A® + (dogl? + dogm® + dyn? 
+ 2dyglm — pv"} AY +- {(dog + das) mn + (dgg + dye) nT} AZ= 0. 
{(dgg + ds) mn + (dys + dz) nl} A” + {(dog + das) mn + (dy + di)nl} 
AY + (deel? + deg? + dg n® + 2dsglm — pv”) A? = 0, 
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The wave equation contains twenty independent constants. 
II. Orthorhombic system (Cay, De, Den). 


(dyP? + dggm® + dzqn* — pv*) A® + (dig + dg) lm AY 
+ (dys + dz) nl A? = 0, 


(dig + dg) lm A® + (dog P+ dygm* + dyyn? — pv) AY 


(dig + dey) nl A® + (dog + dys) mn AY + (degl* + dag” 
+ ds3n* — pv”) A? = 0. 


The wave-equations involve twelve distinct constants. 
Ill. Tetragonal (Cy, S4, Cyn)—tetrad axis parallel to the z-axis. 
{dyl® + dggm® + dygn® + 2d,, lm — pv} A* + {dj (J? — m?) 
+ (dz + dq) lm} AY + {(dsg — dae) mn + (dys + dys) nl} A? =0 
{dy (? — m*) + (diz + deg) Im} A+ {deg ? + dym? + dyn? — 2d,glm 
— pv} AY + {(dis + dys) mn + (dag — yg) nl} A* = 0. 
{(dsg — dyq) mn+(dys + dys) nl} A* + {(dis + dys) mn + (dag — dys) nl} 
AY + { ds; (12 + m*) + dggn? — pv*} A? = 0. 
The number of independent constants contained in the wave equation is fen. 
IV. Tetragonal (Cy, Say, Dy, Dan)—tetrad axis parallel to the z-axis. 


(dy; ? + dggm® + dagn® — pv®) A® + (dig + dgg) lm AY + (dhs 
+ dj) nl A? = 0. 


(dig + dgg) lm A® + (dog ? + dyym* + dygn® — pv*) AY + (dis 
+ d,;) mn A? = 0. 


(dys + das) nl A® + (dys + dys) mn AY + { dj, (J? + m?) + dog n? 
— pv" A? = 0. 


The number of independent constants contained in the wave-equation is 
seven. 


V. Trigonal system (C;, S,)—triad axis parallel to the z-axis. 


{dyl* + dggm* + dg n® + 2d,4mn + 2d,.nl + 2d,4lm — pv} A* 
+ {dyg (m® — I?) — 2d,, mn + 2dygnl + (diz + deg) Im} AY 
+ {dy (7? — m*) + (dag — dae) mn + (dyg+dyg) nl + 2d,5lm}A7=0, 
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{dis (m? — I?) — 2d,, mn + 2d,qnl + (dg + deg) Im} A* 
+ {deg ? + dyym®+ dyn® — 2d,gnn — 2d,, nl — 2d,glm — pv*} AY 
+ {dys (12 — m*) + (dis + das) mn + (dag — gg) nl 

—- 2d lm} A? = 0. 

{dig (2 — m*) + (dyg — dag) mn + (dig + dys) nl + 2dys lm} AX 
+ {dy (2? — m*) + (dis + das) mn + (dag — dye) nl — 2d,glm} AY 
{d,,l* oe dssm? a dygn* -_ pv} A? => 0 

where 
dgg aati (diy ae dis as dso) 


The number of independent constants contained in the wave equations is 
twelve. 


VI. Trigonal (D;, Dsg, Csy)—triad axis parallel to the z-axis. 
{ dyyl? + dgg m? + ds, n® + 2d,4mn — pv*} A® + {2d, nl 
+ (diz + deg) Im}AY + {(dis + das) nl + 2d)s Im} A? = 0. 
{2dy4nl + (diz + deg) lm} A® + {degl? + dym* + dagn* 
— 2d,ymn — pv*} AY + {jg (7? — m?) 
+ (dis + das) mn + dent} A? = 0. 
{(dis + das) nl + 2dys Im} A® + {dys (7? — m?) 
4 (dis + das) mn aa deg nl} AY a“ {d,,/* oe dsm" ok d,3n” 
— pv*} A? = 0. 
where 
deg = (diy — dye — oq) 


The number of independent constants contained in the wave equations is 
eight. 
VII. Hexagonal (Cyn, Ce, Cen)—hexad axis parallel to the z-axis. 
{dy,/* + dss m + dsg n? = 2d. lm — pv"} A*~ + {dis (m? — I?) 
+ (dig + deg) lm} AY + {(dgg — dag) mn + (dis + dys) nl} A? = 0, 
{dig (m® — I?) + (dig + deg) lm} A® + {dog I? + dyym? + dyn? 
— 2d,, lm — pv*} AY + {(dis + dos) mn + (dag — dgg) nl} A? = 0. 
{(dsg — dag) mn + (dys + dys) nl} A® + {(dig + dys) mn 
+ (dig — dg) nl} AY + {dygl? + ds, m? + ds, n® — pv} A? = 0, 


where 


das _ (dy mee dis ee deg) 











64 Sm C. V. RAMAN AND K. S.° VISWANATHAN 


The number of independent constants contained in the wave equations is 
eight. 


VIII. Hexagonal (Dzn, Coy, Dg, Dgn)—hexad axis parallel to the z-axis. 
(d,,)? + dggm* + dss n® — pv®) A® + (diz + deg) lm AY 
ok (dhs + d,s) nl A? = 0. 


(diz + dgg) Im A® + (dggl? + dyym® + dygn? — pv?) AY 
+ (dys + dys) mn A* = 0. 


(diz + dys) (nl A® + mn AY) + (dial? + dg m? + dygn? — pv?) A? = 0, 
Here again the relation 
Age = dy — diz — doo 
characteristic of crystals of the trigonal systems subsists, and number of 
independent constants appearing in the wave equations is six. 
IX. Cubic (T and Typ). 
(dy,l? + dyym® + dssn? — pv?) AX + (dip + dis) (lm AY + nl A?) = 0. 


(dio + dys) Im A® + (dyym? + dyan® + d,,]2 — pv?) AY 
+ (dio —+ ds) mn A? = 0. 


(dip + dys) (nl A® + mn AY) + (dyn? + dal? + dsm? — pv?) A?=0. 


The number of independent elastic constants appearing in the wave equa- 
tions is four. 


X. Cubic (Tg, O and Dp). 


{dyyl® + dy (m® + n2) — pv®) A® + (dy + dys) {lm AY + nl A2) = 0}. 
(dig + dys) Im A* + {dyym® + day (7? + m*) — pv} AY 

+ (dip + ds) mn A? = 0. 
(dio + ds) {nl A*« + mn AY} + {d,,* + dag (? + m?) — pv" A? = 0. 


The number of independent constants appearing in the wave equations is 
three. 


6. STATIC DEFORMATION PROBLEMS 


Equations (1) express the nine stress components in terms of the nine 
strain coefficients. One can work out the inverse transformation of (1) 
and express instead the strain coefficients as linear functions of the stress 
variables. The strain-stress relations can therefore alternatively be written 
also as 


—< , on, (m = 1, 2,....9) (10) 


ne1 
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If 4mn denotes the co-factor of the element dm» in the determinant of 
the transformation (1) and 4 = determinant |dm»|, then Dmn= (4mn/4). 


The strain-energy can be expressed purely as a function of the stress 
coefficients and we have another expression 
2U=7 2 Din TmTn (11) 


for the deformation energy which is equivalent to (3). 


In the classical theory the nine stress variables are reduced to six by 
means of the relations Tyz = Tzy; Tzx= Txz; Txy= Tyx. The deforma- 
tion energy of the solid is thus a function of the six stress variables only. 
Adopting the usual convention of writing T,’, T,’, Ts’, T,’, T;’, T,’ for Txz, 
Tyy» Tzz, Tyz or Tzy, Tzz Or Tzz, Txy Or Tyg respectively, the classical expres- 
sion (4) for the strain-energy becomes 


2U, = ~ a SmnTm’ Tn’ (12) 
The coefficients sm», occurring in the above expression are the well-known 
elastic moduli of the classical theory. One can obtain formule for these 
moduli in terms of our constants Dmy by considering cases of homogeneous 
stresses for which the relations Tyz = Tzy; Tzz = Txz; Txy = Tyzx hold 
good. Making these substitutions in (11) and comparing the resulting expres- 
sion with (2), we get Smn=Dmn when both m and 7 are 1, 2, or 3 and rela- 
tions of the type 


Sia = (Dia + Dys); 

Sag = (Dag + 2D45 + D55); 

Sas = (Dag + Daz + Doe + Doz); (13) 
for the other coefficients. 


We can now write down the expressions for the compressibility, Young’s 
modulus, and the Poisson’s ratio for any crystal in simple terms. In view 
of the fact that all these moduli are determinable from experiments dealing 
with purely homogeneous strains, the formule for them are not essentially 
different from the corresponding ones of the classical theory. By following 
the same methods as those adopted in the latter,* we give below the formule 
for these moduli in our present notation. 


When the crystal is subjected to a uniform hydrostatic pressure P, we 
have from the first three of the equations (10) 
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u, = (Dy + Dy, + D,s) P 
= (Sy) + Sig + S43) P from (13) 

Uz = (Soy + Seq + 59g) P 

Us = (Sg, + Sq + 553) P 
The linear compressibility modulus in the direction of the x-axis is there- 
fore given by 

Ky = (P/uy) = 1/(51. + Sie + 543) (14) 
and the bulk modulus is expressed by the formula 


Kk = Pi(uy + ue + Us) = 1/(Sy. + Seg + S33 + 2523 + 253; + 2549) (15) 
If the solid is subjected to a uniform tension T in the direction of the 
x-axis, obviously T, = T; T,=T;=..T,=0. Hence it follows from 
(10) that the Young’s modulus in the direction of the x-axis is given by 
E = T/u = (1/Dy) = (1/sy) (16) 
Similarly the Poisson’s ratio in the direction of the y-axis is given by 


o = — (Sy2/5y) (17) 
To write down the expressions for the Young’s modulus and Poisson’s 
ratio in any general direction (/, m, n), we require the law of transformation 
of the strain as well as the stress components when one passes over from 
one co-ordinate system to another. If the direction cosines of the axes 
Ox’, Oy’ Oz’ of a new co-ordinate system referred to the original one are 

(1, m, 1), (1, me Mg) and (J; mg ng) respectively, we have 

Tra = 1? Tre + my? Tyy + my? Tzz + myn, (Tyz + Tzy) 
+ ml (Tzx + Taz) + hm, (Tay + Tyzx) 


Tey = hile Tre + mye Tyy + MyM Tzz + (Mmyn, Tyz + man, T2y) 
+ (ml, Tz + Meh Txz) + (hime Tay + lem, Tyx) (18) 


Tya = LbTxx + mm, Tyy + myn, Tzz + (men, Tyz + myn, Tzy) 

+ (Mel, Tze + mle Txz) + (Lm Try + hime Tyx 
Uri? = LPuge + MyWUyy + ny"Uzz + myn, (yz + Uzy) 

+ ml, (Uzx + Uxz) + Lim (Uxy + Uyzx) (19) 
Usty? = Ile Uae + MyM Uyy + MyNe Uzz + (MMe Uyz + MyN, Uzy) 

+ (Mle Uzz + NelUx2) + (Lyme Uxy + ym, Uyx) 


Uy? = Ile Ua + MyMg Uyy + MyNg Uzz + (Mey Uyz + MyNg Uzy) 
+ (Ml; Uae + Mylguaz) + (lnywuay + Lmeuyz) 
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Now, if the solid is subjected to a uniform tension T over the planes x’ 
= const. of which the normal has direction cosines (/,m,,), we get from 
the transformation rules 


Txx = LT: Tyy,= m?2T; Tz = n,*T; Tyz = Tzy => mynyT ; 
Tze = Taz = mhT; Try = Tyr = lym,T 
Further, uzz is given by (19). Hence denoting by qy the row vector 


(1? m,? n? m,n, nl, 1m,) with six components, we find the expression 
for the Young’s modulus in the direction of the vector (/, m,7,) as 


Ey = 1/qx S qx’ (20) 


where S denotes the matrix (smn). Similarly the Poisson’s ratio in the direc- 
tion of the y’-axis which has direction cosines (/, mz n2) with respect to Oz, 
O, and Oz is given by 

eae ag 2 6 


2 o¢ op ap 


0 = — 3g [Mae 1 ame 1 ans + Gm) +4 SEHD 


+ him <7 (21) 


where ¢ = qy SGy and the differential coefficients are formed as if these 
arguments are independent. 


It is interesting to note that only twenty-one of the forty-five constants 
figuring in (10) appear in the above formule. As mentioned earlier, this 


is a consequence of the homogeneity of the strains applied to evaluate these 
static moduli. 


The stresses and strains which appear when elastic materials are subject 
to torsion or flexure are essentially heterogeneous, and hence in dealing with 
them, our elastic moduli Dm, will appear in combinations other than those 
which figure in homogeneous deformations. The present theory is quite 
competent to handle such problems, but to deal with them in detail would 
carry us far beyond the scope of this paper. 


7. SomMe CONCLUDING REMARKS 


The main purpose of the present paper has been to establish the neces- 
sity for an amendment of the phenomenological theory of elasticity as 
universally accepted hitherto. Any theory to be acceptable should include 
in its scope elastic stresses and strains of the most general type and specify 
them in an analytical form from which the equations of equilibrium in static 
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problems and the equations of motion in wave-propagation can be written 
down immediately. These requirements are not met by the present form 
of the theory but are completely satisfied if all the nine components of stress 
and strain in the usual tensor formulation are retained. When thus amended, 
the phenomenological theory is capable of handling all the problems of the 
subject, including those which arise in its practical applications. The more 
important of these are considered in detail in the course of the paper and the 
results are set out explicitly for the different crystal classes so that they could 
be readily made use of. 


As already remarked in the introduction, the phenomenological theory 
as amended gives results in complete accord with those derived from the 
atomistic approach to the theory of elasticity of crystals based on the most 
general scheme of interatomic forces. To discuss the latter further or to 
give an account of the somewhat confused history of the subject would lie 
outside the scope of the present paper. A few remarks regarding these 
matters will however be found in the attached appendix. 


SUMMARY 


The fundamental aspects of the phenomenological theory of elasticity 
are critically examined and it is shown that the tensor representation of the 
elastic strains and stresses in the general case should be in the unsymmetrical 


form. On this basis, the stress-strain relationships are deduced and tabu- 
lated for the different crystal classes. The equations determining the veloci- 
ties of wave-propagation in different directions are also obtained and tabu- 
lated. Static deformation problems are then discussed and it is shown that 
in the particular case of homogeneous strains, the elastic constants group 
themselves in linear combinations which are equivalent to the elastic modulii 
of the theory in its familiar form. In wave-propagation, however, the strains 
and stresses are heterogeneous and hence all the elastic constants are involved 
and appear in linear combinations which are different and also larger in 
number than those which figure in the formule for homogeneous deforma- 
tions. These results are completely in accord with the consequences of 
the atomistic theory based on interatomic forces of the most general type. 
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APPENDIX 


As is well known, the theory of the elasticity of solids in its present 
form was initiated by Cauchy who put forth the fundamental idea of expressing 
elastic stresses and strains in the manner adopted in the phenomenological 
theory, and also developed the theory on an atomistic basis. Later theorists, 
notably Green, Stokes and Kelvin retained Cauchy’s method of specifying 
elastic stresses and strains but adopted a purely phenomenological approach. 
The view advocated by them that a completely zlotropic body would have 
twenty-one elastic constants and not fifteen as derived by Cauchy received 
general acceptance. 


In a paper published in these Proceedings some years ago by one of 
us (Raman, 1943), a theory of the dynamic behaviour of crystal lattices was 
developed based on the assumption of interatomic forces of the most 
general type. In subsequent years, the consequences of that theory was 
worked out in detail for the case of diamond and confirmed by a series of 
spectroscopic investigations on the scattering of light, the luminescence and 
infra-red absorption by that crystal. The high values of the interatomic 
force-constants disclosed by the spectroscopic behaviour of diamond were 
evidently related to its exceptional elastic behaviour. In the endeavour to 
place this relationship on a quantitative basis, a fresh approach was made 
by one of us (Viswanathan, 1954) to the atomistic theory of the elasticity 
of crystals, both from the static and dynamic points of view, and some sur- 
prising results emerged. Contrary to the assertion made in the papers of 
Max Born and his collaborators which also finds a place in their recent book 
(Born and Kun Huang, 1954), no difficulty was encountered in expressing 
the energy of static deformations in terms of interatomic forces of the most 
general type. It was found that this expression contained forty-five inde- 
pendent constants, but for homogeneous or irrotational strains they appeared 
in twenty-one distinct linear combinations. The dynamics of wave-propaga- 
tion in crystals was also investigated and it was shown that the expressions 
for the wave-velocity contained the same forty-five constants but in different 
linear combinations. The work of Born and his school on the dynamic 
problem was critically examined and it was shown that the assumptions 
made by them in the attempt to reduce the forty-five constants which appeared 
in their theory to twenty-one had no theoretical justification. 


The results of Born and his school were also contradicted by Laval in 
some recent publications (1951). More recently still, a series of papers 
have been published by Le Corre in which Laval’s ideas have been further 
developed, On reading those papers, one obtains the impression that their 
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author believes the results of the atomistic and phenomenological approaches 
to differ essentially. For example, in order to account for the non-sym- 
metric character of the stress tensor, internal couples are postulated to 
balance the differences in the angular momenta. We may remark that there 
is no room for such a postulate, since the analytical specification of the stresses 
in terms of the tensor components should itself suffice to describe the state 
of the solid completely. In the last paper of the series, a statement also 
appears that an atomistic approach is essential to solve such familiar problems 
in elasticity as torsion and flexure. 


In conclusion, we have to thank Mr. A. K. Ramdas for his help in 
the preparation of the tables appearing in the paper. 
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ABSTRACT 


As a part of an attempt to develop new methods of measurement 
of diffusion constant of metals and non-metals into minerals using 
radioactivity techniques, expressions are derived to relate the processes 
of diffusion and the emission of alpha-particles under three sets of possible 
experimental conditions. The first case pertains to a plate whose thick- 
ness is equal to the characteristic range of the alpha-rays of the radio- 
nuclide used. Under the second case is considered a plate thicker than 
the characteristic range of the alpha-rays in the plate. The third deals 
with cases wherein the radioactive matter used emits alpha-rays with 
different ranges less than the thickness of the plate. 


INTRODUCTION 


THE problem of diffusion of metals in the solid state has been extensively 
studied during the last three decades (Mehl, 1936), because of its importance 
in the industrial operations like the preparation of special alloys and the 
development of nitriding, chromizing and calorizing, sherardizing and silico- 
nizing and the formation of bi-metal strip and veneering metals. Diffusion 
data are also essential for studies on crystal physics and the elucidation of 
phenomena such as annealing, age hardening, plasticity and order-disorder 
transformations in alloys. Our present state of knowledge on the subject 
of solid diffusion has recently been reviewed by Barrer (1951). 


The voluminous literature on metal diffusion stands in marked contrast 
to the virtual absence of any systematic work on the measurement of diffu- 
sion constants of metals and non-metals in minerals and rocks. The dis- 
covery of radioactive isotopes gave a great fillip to self-diffusion studies on 
metals (Hevesy et al., 1913, 1925, 1929 a, b, 1932, 1933; Groh and Hevesy, 
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1920; McKay, 1938; Steigman etal., 1939). The pioneering work of 
Jensen (1952) on the solid diffusion of radioactive sodium (Na??) in perthite 
shows tremendous potentialities for further researches on the subject. 


The principles of diffusion have been discussed in detail in an earlier 
communication (Aswathanarayana, 1954). As a part of the long-range 
programme of theoretical and laboratory investigations on the measure- 
ment of diffusion constants, an attempt is made now to derive expressions 
relating the processes of diffusion and emission of alpha-particles. 


GENERAL STATEMENT OF THE PROBLEM 


Let ABC D-A’ B’C’ D’ with dimensions a, b, c (a being the thickness) 
represent a thin, uniform section of a mineral. A uniform, thin coating 
of an alpha-emitting radio-nuclide is given to the face ABCD. Let v be 
the concentration of radioactive matter on that surface. Simultaneously 
with the above operation, an a-particle detecting device—nuclear emulsion 
plate suitable for alpha-track radiography—is kept in contact with the face 
A’ B’C’ D’ of the thin section. Each radioactive atom undergoing decay 
acts as a centre of a sphere of radius equal to the range of alpha radiation 
in the medium. When the sphere intersects the plane of the emulsion, it 
leaves behind its imprint in the form of an image (Yagoda, 1949). The 
number of alpha-particles emerging from the lower face A’ B’C’ D’ and 
getting recorded on the nuclear emulsion plate is a function of (1) the con- 
centration (v) of radioactive matter on the face ABCD, which determines 
the rate of production of alpha-particles, (2) the characteristic range (p) 
of the alpha-particles of the radioactive nuclide in the mineral and (3) the 
thickness of the plate (a). With time, a part of the material of the radio- 
active nuclide diffuses into the mineral and this results in a larger number 
of alpha-rays emerging from the face A’ B’C’D’ of the mineral section 
and getting recorded on the plate. Thus from a knowledge of (1) the range 
(p) of alpha-particles of the radioactive nuclide, (2) the concentration (v) 
of radioactive matter on the plate and (3) the number of alpha-particles 
(N) recorded on the plate after a time T, it is possible to determine the diffu- 
sion constant (D). 

MATHEMATICAL TREATMENT 
The following assumptions are made in this context:— 


(i) The diffusion of radioactive matter in a mineral is governed by laws 
analogous to laws of conduction of heat in a solid and (ii) the mineral exhibits 
diffusion isotropy. 
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The concentration v (x, y, z, ft) is governed by 


ov 2 y 
= DV 


The experimental conditions may be taken as equivalent to the diffusion 
of matter in a rectangular slab of thicknesses a, b, c, where ais the distance 
between the plane on which radioactive matter has been spread and the 
opposite plane (which is in contact with the nuclear emulsion plate). 


The boundary conditions are: 

at x=0 v=, 
x=@a 
y= 0 
y=6b 
z=0 
z= Cc 

at t=0 v= 0. 

Then v can be written as the sum of two terms U and W (say)—U 


being a function of x, y, z alone and W, a function of x, y, z, t (Carslaw and 
Jaeger, 1947, pp. 18-19, 165). 


Concentration v = v (x, y, z, t) 
6 co sinh a (a — x) sin (2p ah ™Y sin =a 
a 
v (x,y, 2,0) = > ya (2p + 1) (2q + 1) sinh aa 


lax mry nnz 


Je-Po", m, n, t e aie 
_ 3 >I >) - le sin — = sin b sin = 
mn? mn? nea? 
1=1 


wr'-+ a “—_ me \ sinh a = aie 








where 
_ p+ Na, Cg+ Ia 
“pb? 


[2 m? n2 
— bo 
a7) m= 7 (2 + b2 + cz) ° 


Case I.—The following experimental conditions are assumed in this 
case: (i) the thickness (a) of the slab is equal to the range (p) of the a- 
particles of the radio-nuclide, (ii) the a-rays are mono-energetic or in other 
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words, they all have a characteristic specific range (p) in the mineral, 
(iii) the radioactive matter has been spread on the mineral surface uniformly 
and (iv) the nuclear emulsion plate is in contact with the lower side of the 


mineral plate simultaneously with the spreading of radioactive matter on 
the upper surface. 


The number of a-particles emitted between the time instants 0 and 
T is given by 


N = fat f dz f dy fords 


where A is the fraction of radioactive matter that disintegrates per 
second. 





Putting 
_ 64 ,be 
A => ma, ‘sinh (aya) [cosh (aga) <a 1], 
a 256 v,ab>c* : 
n® {b2c? + a? (b? + c%)} sinh fe JP +c 
a; 
aj —7 ale + 2 ’ 
- . =. eS 
e ~eal> + Bt 
we get 


2T 2 
— = (A — B) T+ — D (approximately), 


taking into account the terms in v with p = gq = 0, and ] = m =n = 1, 
From this equation D can be calculated. 


Case II.—In the second case, the following experimental conditions 
are assumed: (i) The thickness (a) of the plate is larger than the range 
(p) of the alpha-particles of the radio-nuclide in the mineral, (ii) the alpha- 
particles are mono-energetic, (iii) spreading of the radioactive matter on the 
mineral face A BC D is uniform and (iv) the nuclear emulsion plate is in 
position by the time radioactive matter reached x = a—p plane. 


T c b a 
> =f at f dz f dy { vdx 
5 t) 0 a-p 











af 
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where 7 is the time taken by the readioactive matter to reach the plane 
x=a-—p. It is given by 


;' mil + a (Fa + a) \ sinh 7a vp +2 sib (agp) 
i at ae 
2 sinh (aoa) sin 7 (1 ) 





approximately. 


Putting 


64v,bc 


A’ = 7 sinh (aya) [cosh (aop) — 1], 


p = : {1 + cos7(1 -°\ a? (E* — 77), 


y=A'—3 {1 + cos «(1 ~= 


(B, a, a) the same as above in Case I), we get 


“ = y(T—7)+ 6D (approximately), 


whence D can be calculated. 


Case III.—This case is similar to the second case but for the assump- 
tion that the alpha-particles emitted by the radioactive matter are not mono- 
energetic, i.e., they have got several ranges as different from a single range 
considered under the second case. 


Let p, p’, p’....p\™ be the ranges of a-particles emitted by the radio- 
active matter. Let p <p’ < p’....p\™ <a. Let 7” be the time taken 
by the radioactive matter to penetrate to such depth that radiations from 
it of range p'” can influence the photographic plate. Let T be the time 
that has elapsed upto the removal of the photographic plate. 


Then 
Eg e b a 
Da fa fae fay f vas 
T 0 0 a—p 
Ky a 


+ fat ff de [% J vdx, 


nin) a—pm 
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p=p* and r= 7"), 


Nu 
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Let the values of 8 and y in Case II be termed as B“ and y‘*’ when 


N 


Then 


_= 2 ys(T — 7s) + D 2 Bs (approximately). 


From this formula D can be calculated. 
Further investigations are proceeding. 
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1. INTRODUCTION 


NickEL sulphate crystallizes from solution at room temperature (30° C.) 
as NiSO,.7H,O in the orthorhombic system. At about 35° C it crystallizes 
as NiSO,.6H,O in the tetragonal system, while at about 60° C. it crystallizes 
as NiSO,-6H,O but in the monoclinic system. The complete Raman 
spectrum of crystalline nickel sulphate has not so far been reported since 
the 44358 and A 4046 radiations of the mercury arc ordinarily employed 
are absorbed. Nisi (1930) using the group of lines at 43126 in the mercury 
arc spectrum has reported in the case of NiSO,.7H,O the principal sulphate 
frequency of shift 985cm.-! and two water bands at 3209 and 3432 cm." 


The absorption in the near ultraviolet for nickel sulphate lies at A 3850 
and the crystals are transparent to radiations of shorter wavelengths over 
a wide range in the ultraviolet. The present paper reports the results of an 
investigation with the orthorhombic and tetragonal crystals using the reso- 
nance radiation A 2536-5 of the mercury arc for exciting the Raman effect. 


2. EXPERIMENTAL DETAILS 


The crystals of NiSO,.7H,O and NiSO,.6H,O were grown by the 
method of slow evaporation, the latter at a temperature of 35° C. in a thermo- 
stat. The crystals of NiSO,.6H,O were of a deeper green in colour than 
those of NiSO,.7H,O. The biaxial nature of the interference figures in 
the case of NiSO,.7H,O and the uniaxial nature of the figures of NiSO,.6H,O 
observed under the polarising microscope confirmed their respective identi- 
fication from depth of colour and external form. The crystals of 
NiSO,.7H,O were in the form of rectangular tablets of size 20 mm. x 10 mm. 
x2mm. The crystals of NiSO,.6H,O which were pyramidal in form were 
rather small. Since both the crystals were efflorescent, a thin film of glycerine 
was put over thera to reduce the rate of efflorescence. However, the crystals 
of NiSO,.7H,O deteriorated rapidly and fresh crystals had to be used fre- 
quently. 


The Raman spectra were recorded with a Hilger medium quartz spectro- 
graph having a dispersion of 140cm.-' per mm. in the \2536-5 region. A 
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water-cooled magnet controlled quartz mercury arc served as an intense 
source of A 2536-5. The crystals were placed near the arc and the scattered 
radiations were focussed on the slit of the spectrograph with a quartz con- 
denser. The scattered resonance radiation was filtered out by the mercury 
vapour inside the spectrograph. With a slit width of 0-04-0-05 mm.,, 
exposures of the order of ten hours were given to obtain reasonably intense 
spectrograms. 
3. RESULTS 


Figs. (1) and (2) in the accompanying Plate reproduce the spectra 
obtained with NiSO,.7H,O and NiSO,.6H,O respectively. The lowermost 
part of the spectrum in Fig. 2 which is very intense due to parasitic light 
helps to identify the lines of the mercury arc. 


An examination of the spectra of the two crystals reveals the following 
features: (1) The low frequency spectra of these crystals are faint and 


TABLE I 
Frequency Shifts in cm. 





NiSO,-7H,O NiSO,.6H,O 





~ 75 v.w. 
~ 110 v.w. 
Lattice ~ 156 w. 


Oscillations 210 m. 210 m. 
255 b.m. 236-69 b.m. 


400 w. 
444 s. 440 s. 
465 s. 469 s. 








620 ? 620 ? 





986 v.s. 988 v.s. 











1060 b.s. 
1098 b.m. 1092 b.s. 
% 1138 b.s. 1133 b.s.} 
~1165 v.w. 








Water Bands ~3282 v.s. ~3257 v.s. | 
| ~3438 v.s. |~3439 v.s. 


| 








| 

s., Strong; w., weak; v.S., very strong; v.w., very weak; m., medium; b., broad. 

~ value of the frequency shift estimated from the relative position of the line with respect to 
the neighbouring iron arc lines, 
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diffuse. (2) The frequency shifts due to the oscillations of the sulphate 
ions in the two crystals are different in the two cases. It may be noticed 
that the splitting of the doubly degenerate oscillation at 450 cm.“ is larger 
in the case of NiSO,.6H,O than in the case of NiSO,.7H,O. Measure- 
ments indicate that the non-degenerate oscillation has a frequency shift of 
986cm.! in the case of NiSO,.7H,O, while in the case of NiSO,.6H,O 
it has a shift of 988 cm. It will also be noticed that the triply degenerate 
oscillation at 1100cm.-' has split into more components in the case of 
NiSO,.7H,O which is of lower symmetry than NiSO,.6H,O. The water 
bands also show differences in their structure in the two cases. 


Table I gives the frequency shifts observed in the two cases. The fre- 
quency shifts of some of the lines which were faint and diffuse could only 
be approximately estimated. The sulphate ion frequency of shift about 
620cm.-' could not be unambiguously observed owing to the proximity of 
the intense mercury line at A 2576. 


4. DISCUSSION 


NiSO,.7H,O has been assigned to the space group P 2,2,2, (D,*) of the 
enantiomorphous hemihedral class of the orthorhombic system, whereas 
NiSO,.6H,O has been assigned to the space group P 4,2,2(D,‘) of the 
enantiomorphous hemihedral class of the tetragonal system. (Beevers and 
Schwartz, 1935; Beevers and Lipson, 1932). In the case of NiSO,.6H,O 
the nickel and sulphur atoms occupy the special positions with the symmetry 
of the two-fold axis of rotation. The six water molecules are said to be 
approximately octahedrally co-ordinated around the nickel ions. On the 
other hand, in NiSO,.7H,O the nickel and sulphur ions do not occupy any 
such symmetric special positions. Six of the seven water molecules are 
approximately octahedrally co-ordinated around the nickel ions in this case 
as well. In the case of NiSO,-7H,O as well as in NiSO,.6H,O there are 
four molecules per unit cell. 


The degradation of the symmetry of the sulphate ion in the crystalline 
state and the consequent removal of the degeneracy of the oscillations results 
in the splitting of the degenerate oscillations, and hence all the nine frequencies 
of the SO, ion should appear. The coupling between the four molecules 
in the unit cell further increases the number of these internal oscillations. 
Thus, in the case of NiSO,.7H,O thirty-six components should appear 
whereas in the case of NiSO,.6H,O this number is only twenty-four owing 
to the higher symmetry of the crystal. Though the theoretically expected 
multiplicity of lines is not fully manifested in the spectra, the differences in 
the nature of the spectra due to the difference in the crysial symmetry and 
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tructure are evident. The breadth of some of the lines (e.g., the group near 
1100 cm.-') might be attributed to the overlapping of several of the compo- 
nents and their inadequate resolution. 


The Raman spectrum of NiSO,.7H,O exhibits a close correspondence 
with the spectrum of MgSO,.7H,O. (Shantakumari, 1953) the crystal struc- 
ture of which is isomorphous with that of NiSO,-7H,O. The following 
frequencies have been reported by her. 56, 75, 97, 118, 150, 252, 303, 445, 
461, 609, 620, 986, 1059, 1064, 1076, 1098, 1134, 1148, 3185, 3228, 3338, 
3406, 3446 and 3485 cm.-! 


Schaefer and Schubert (1916) report from their studies on the infra-red 
reflection spectrum of NiSO,.6H,O,. frequencies at 635, L111, 1116, 2940, 
3030 and 3125cm.-! No infrared absorption studies appear to have been 
made with crystalline nickel sulphate. 


In conclusion, the author wishes to express his sincere thanks to 
Prof. Sir C. V. Raman, F.R.S., N.L., for the kind interest that he took during 
the course of this investigation. 


5. SUMMARY 


The Raman spectra of orthorhombic NiSO,.7H,O and the tetragonal 
NiSO,.6H,O have been investigated using the resonance radiation of mer- 
cury for excitation. Readily observable differences are exhibited by the 
spectra in the two cases which may be ascribed to the differences in crystal 
structure and symmetry. 
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1. INTRODUCTION 


A stupy of the thermal expansion of calcite over a wide temperature range 
is of interest for the following reasons. Calcite belongs to the trigonal 
system and hence possesses two principal expansion coefficients—one along 
the optic axis and the other at right angles toit. Also its anisotropy of thermal 
expansion is large as it is the general case for crystals of the type M (XO,) 
(Lonappan, 1955). Its expansion coefficient perpendicular to the axis is 
negative and small. Such a negative expansion coefficient is of rare occur- 
rence. Also its Raman and infrared spectra are known in great detail 
(1945, 1953). The temperature variation of its Raman spectrum has been 
investigated by Narayanaswamy (1946). Bhagavantam and Venkatarayudu 
(1939) have calculated the frequencies of all the normal modes of the crystal 
lattice employing a valence force field. Yet it is definitely not known whether 
calcite is typically ionic or typically of the valence type. X-ray studies 
indicate the former type; while its intense Raman spectrum is what is to be 
expected from a crystal belonging to the latter type. Chapman, Topping 
and Morral (1926) and Lennard Jones and Dent (1927) assumed a hetero- 
polar force field and derived its lattice constant and other properties. An 
attempt at explaining the thermal expansion data in terms of the Raman 
and infrared spectra will be very interesting. The elastic constants of calcite 
have been measured by Bhimasenachar (1945). The specific heat, dielectric 
constant and a host of other physical properties have all been determined 
for calcite. 
2. PREVIOUS WORK 

The thermal expansion of calcite has been measured near room tempe- 
rature by Fizeau (1868) and Benoit (quoted in Weigle and Saini’s paper) 
using the interference method. The expansion of calcite down to liquid 
hydrogen temperatures was measured by Adenstadt (1936) using the inter- 
ference method. 


The lattice expansion of calcite was measured by Megaw (1933) from 
room temperature to 100°C. She claims an accuracy of 10% and her 
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values are in good agreement with the macroscopic expansion values of 
Fizeau. Weigle and Saini (1934) determined the lattice constant and thermal 
expansion of calcite using its powder photograph taken with a Seeman- 
Bohlin camera. They used a crystal which was not very perfect as they felt 
that thermal expansion might not be affected by small impurities or lattice 
imperfections. They claim an accuracy of better than 1%. Yet their values 
for expansion between 0 to 100°C. are lower than those of Fizeau and 
Benoit. The discrepancy is large—about 20%. Yet they conclude that 
this should not be taken as evidence for any real difference between lattice 
and macroscopic expansion. Such a discrepancy might perhaps have been 
caused by the imperfections in the crystal. Recent experiments by Glover 
(1954) on the lattice expansion of pure and calcium-containing potassium 
chloride have definitely shown the influence of impurities to be negligible. 
In view of this finding, Weigle and Saini’s data appear to be unexplainable. 


No measurements of high temperature expansion are present and the 
present investigation was carried out to fill this lacuna. 


3. EXPERIMENTAL DETAILS 


The experimental arrangement employed was that set up by Press (1949). 
The method of measurement has already been described in a previous paper 
by the author (1955). Preliminary measurements of the expansion coeffi- 
cients both parallel and perpendicular to the axis were made by Mr. A. K. 
Sridhar; but he obtained anomalous results, probably because the orienta- 
tions of his specimens were incorrect. The same specimens were used in the 
present investigations the orientation having been corrected with a polarising 
microscope. The specimens used finally were correct to within a degree of 
the required orientation. Three specimens were cut with the optic axis 
parallel to their bases. These were used to measure a,. Three others were 
cut with the optic axis perpendicular to their bases. These served to deter- 
mine a,. 

4. RESULTS AND DISCUSSION 
The results of these measurements are given in Table I below. The 


curves showing the variation of a, and a, with temperature are given in 
Fig. 1. 


The pieces employed for determining a, were only 0-312 cm. thick and 
as the expansion coefficient is very small, reliable measurements for the 
mean expansion coefficient over as much as 100° C. only could be obtained. 
Even then the accuracy of measurement is reduced and is about 5% in this 
case, 
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TABLE I 
Values of a, and a, for calcite at various temperatures 


























Range of Mean Range of Mean 
Temp. Temp. a4, X 10° Temp. Temp. a. X 10° 
AT°C. oe rk ng On eG 
64-85 74:5 27-2 45-115 80 —5-78 
85-105 94-9 28-2 70-150 100 —5-00 
124-143 133-4 29-2 99-187 143 —4-57 
179-198 188-9 29-9 143-241 192 —4:-09 
217-235 225-6 31-2 187-289 238 —4-09 
250-268 259-3 31-4 241-347 294 —4-04 
296-314 304-9 32°4 310-420 365 —3-83 
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The curves for a, and a, as determined by the author are shown in 
Fig. 1. The low temperature expansion values of Adenstadt (1936) are 
plotted as crosses in the same figure. It is seen that the low temperature 
values of Adenstadt and the high temperature values of the author lie on 
a single curve for both a, and a,. 


The values of a, and a, as determined from the graph at 50°C. are 
compared with the values in literature. 


a, at 50°C. = 26-6x 10 
a, at 50°C. = — 52x10 


The values in the literature are 


a, = 26:2x 10° a, = 25-3+1-5x10°8 
a, = — 5-40%10-8} Fizeau q) = — 5-2 £O-7x10° Megaw 


a, = 25-7x 10-8 21:0x10°* | 


. a . es 
oe 5-50 10-* Benoit a, = — 3-8x 10-6 j Weigle and Saini 


The agreement with the values of Fizeau, Benoit and Megaw is very 
good. The values of Weigle and Saini are too low presumably because of 
the imperfect crystal they employed. 


In the following table are collected the thermal expansion values a, 
and a, at various temperatures, the volume expansion coefficient B, the speci- 
fic heat per gram cp and the Gruneisen constant y. The values of cp have 
been taken from International Critical Tables; the value of the compressi- 
bility X is taken as 15-4 10-!* sq.cm./dyne from Bhimasenachar (1945). 

TABLE II 
The linear expansion coefficients a, and a,, the volume expansion coefficient B, 
the specific heat and the Gruneisen constants for calcite at 


different temperatures 
Compressibility X =1-54x 10-2”. Molar Volume V =34-2 c.c. 





Se Gruneisen 
rc a4, X 10® a,x 108 Bx10®  joules/gm. Const. y 





16-4 —6°8 2°8 0-452 0-14 
20-6 —6°3 8-0 0-590 0-30 
25-0 —5°4 14-2 0-763 0-41 
28-1 —4-9 18-3 0-870 0-47 
30°8 —4:2 22°4 0-962 0-52 
32°8 —3-9 25-0 1-046 0-53 
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As usual the Gruneisen constant decreases with decreasing temperatures. 
The small value of the Grunesien constant is noteworthy as it is typical for 
a valence crystal. The intense Raman spectrum of calcite and its dielectric 
polarisation also point to such a type of binding in the crystal. 


The author is indebted to Prof. R. S. Krishnan who suggested this 
problem and provided the specimens. He is also thankful to Dr. Vedam 
for his help in determining the orientation of the specimens. 


5. SUMMARY 


The thermal expansion of calcite along and perpendicular to the optic 
axis were measured from room temperature up to 400°C. a, increased 
with increasing temperature, while a,, which is negative, decreases in magni- 
tude with increasing temperatures. The high temperature curves join up 
well with the low temperature curves of Adenstadt. The Gruneisen constant 
is shown to vary with temperature in the usual manner. Its small value 
points to a valence type of binding in the crystal—a fact supported by the 
intense Raman effect and dielectric polarisability of calcite. 
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1. INTRODUCTION 


ABSORBING biaxial crystals in general display a variety of remarkable optical 
phenomena in the vicinity of both the optic axes. For example, if an 
extended source of unpolarised light is viewed through a plate of highly 
pleochroic material cut normal to an optic axis, two dark brushes—the 
Brewster’s brushes—are generally seen in the field of view; while if a polariser 
be inserted in front of the plate, the so-called idiophanic rings are observed— 
similar to the interference rings that can appear in the case of a transparent 
crystal if an analyser be also present. 


The theoretical investigations of Waldemar Voigt focussed attention on 
the fact that certain of the features relating to the propagation of light in 
the vicinity of an optic axis differ radically from those obtaining in trans- 
parent media. Thus, whereas along any general direction in a transparent 
crystal there are two particular linearly polarised vibrations that can be 
propagated without change of form, this is no longer the case in absorbing 
crystals. As a matter of fact, close to an optic axis and on either side of it, 
there even exist two directions—the singular axes—with the following re- 
markable properties: only a right-circular vibration can be propagated 
without change of form along one of these axes, and only a left-circular 
vibration along the other.* In this paper it will be shown that the various 
features of the propagation of light in absorbing media may also be con- 
veniently regarded as due to the superposed effects of birefringence and 
dichroism. Because of the simplicity of the method, it has also been possible 
to make a more detailed investigation of the following interesting question: 
what will happen when, for example, a right-circular vibration is incident 
in the direction of a singular axis where only a left-circular vibration can be 
propagated without change of form? The results obtained in this connection 
are at variance with those expected by Voigt (Section 6). 


*A non-mathematical summary, in English, of the main results of Voigt’s investigations 
may be found in Reference 1. 
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The analysis of the propagation of light in absorbing biaxial media from 
the standpoint of the electromagnetic theory has been discussed by several 
authors. The comparatively simple case of orthorhombic crystals (where 
the principal axes of the dielectric tensor and the complex conductivity tensor 
necessarily coincide) is discussed in Drude’s treatise.? But the somewhat 
oversimplified presentation there given omits entirely those theoretical and 
experimental features with which we shall be particularly concerned. These 
may be found described in more detailed treatments,*~* particularly those 
of Voigt and Pockels. In view of the complexity of the phenomena involved, 
it would appear that a consideration of the problem from a simpler though 
less rigorous standpoint would certainly be useful. Such an approach is 
provided in the present paper, and, as we shall show, the method adopted 
leads to results that are practically identical with those of the electromagnetic 
theory. 


2. GENERAL FEATURES OF LIGHT PROPAGATION IN ABSORBING MEDIA 
(a) The Index and Absorption Ellipsoids 


In an absorbing biaxial medium not possessing optical activity, the two 
waves propagated along any direction appreciably inclined to both the optic 
axes may be regarded as practically plane polarised (though not rigorously 
so, as in the transparent crystal). And as in a transparent crystal their 
‘ vibration-directions ’ may then be considered to lie on the principal planes, 
their velocities being determined by their vibration-directions thus: the 
reciprocal of any radius of a so-called index ellipsoid gives the velocity for 
vibrations parallel to that radius. In addition, the two waves have different 
coefficients of extinction «, and «,, these being determined again by their 
vibration-directions thus: the reciprocal of any radius of a so-called absorp- 
tion ellipsoid gives the value of 4/(2«v*/c) for vibrations parallel to that 
radius, v being the velocity for that vibration-direction. 


By assuming that the above statements hold good even for directions 
in the vicinity of an optic axis, it is indeed possible to explain some of the 
phenomena observed there—and such a procedure is in fact followed in 
Drude’s treatise. For example, the occurrence of Brewster’s brushes can 
be explained along the following lines. In the neighbourhood of an optic 
axis, a comparatively small change in the direction of propagation will in- 
general cause an appreciable change in the inclinations of the two principal 
planes to the axial plane; this in turn will lead to a large variation in the 
total absorption, since the absorption coefficients of the two waves will be 


determined by the orientation of their vibration-directions. 
A3 
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(6) The Elliptical Polarisation of the Waves 


The appearance of idiophanic interference rings with a polarizer alone 
cannot however be explained on the assumption that the light incident along 
any direction is split up into two linearly polarized beams with their vibra- 
tions at right angles to one another; their states of polarization being ortho- 
gonal, two such beams will be incapable of interference with one another 
(unless brought to the same plane of vibration by an analyser). 


The fact is that when we turn to directions of propagation in the vicinity 
of an optic axis, we are no longer justified in neglecting a remarkable and 
important consequence of the phenomenological theory: namely, that the 
two waves propagated in any general direction in an absorbing biaxial medium 
are in reality, elliptically polarized. Though the two elliptic vibrations have 
their axes majores at right angles, and their ellipticities equal, they are 
rendered non-orthogonal by the fact that they are of the same handedness; 
and this last mentioned feature (together with the fact that the major axes 
do not in general coincide with the principal planes) distinguishes the situa- 
tion sharply from that obtaining in optically active (transparent) crystals. 


In the context of the elliptical polarisation of the waves, the index and 
absorption ellipsoids—strictly speaking—retain significance only in terms 
of the dielectric and conductivity-like tensors by means of which they are 
defined. Nevertheless, as we shall show, the existence of the two non- 
orthogonal elliptically polarised waves may be conveniently treated as due 
to the superposed effects of birefringence and dichroism—just as the propa- 
gation of two orthogonal elliptically polarised waves near to an optic axis 
in an optically active transparent medium, may (by Gouy’s hypothesis**) 
be conveniently treated as due to the superposed effects of birefringence 
and rotation. 


3. THE SUPERPOSITION OF BIREFRINGENCE AND DICHROISM 


Consider a plate cut perpendicular to an arbitrary direction z which is 
also taken as being normal to the plane of the paper. Let OX; and OY; 
(Fig. 1) be the trace of the principal planes of refraction—defined as usual, 
either in terms of the index ellipsoid or the optic binormals. Similarly let 
OX;, and OY;, be the trace of the principal planes of absorption—which we 
shall define analogously, either as containing the major and minor diameters 
of the elliptical section of the absorption ellipsoid made by the plane of the 
paper, or as the internal and external bisectors of the angle subtended on 
the z-direction by the two absorption-binormals (normals to the circular 
sections of the absorption ellipsoid). 
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P, and P,—States of polarisation propagated unchanged. 
OX;,, and OY;—Trace of the principal planes of absorption. 
OX, and OY,—Trace of the principal planes of refraction. 


Consider an arbitrarv’ elliptic vibration P (which, for the sake of con- 
creteness, may be temporarily identified with the one marked Py in the figure), 
In the absence of absorption the arbitrary elliptic vibration P will be resolved 
into two vibrations (along the principal planes of refraction, OX; and OY;) 
between which an infinitesimal phase difference dz will be introduced corres- 
ponding to a passage dz. Since, however, anisotropic absorption is also 
present, we perform, in addition, the infinitesimal operation of linear 
dichroism; the elliptic vibration—as modified by the infinitesimal operation 
of birefringence—is resolved into two linear vibrations (this time, along the 
principal planes of absorption, OX; and OY ,), the amplitudes of which 
are then reduced by the multiplying factors (1 — k,dz) and (1 — k,dz) res- 
pectively. The differential absorption of the two components will cause 
the state of the elliptic vibration to ‘ move towards’ the state of polarisation 
of the less absorbed component OX;, (a phrase which acquires a more vivid 
meaning in the Poincaré sphere representation). Those states of polarisa- 
tion alone can be propagated without change of form, which under the 
successive infinitesimal operations of birefringence and dichroism (applied 
in either order) remain unaltered in form and orientation—and to these 
states of polarisation alone can definite velocities of propagation and 
coefficients of absorption be assigned. 


Several particular cases may first be noted. Along the optic axial 
directions where the birefringence vanishes, the two waves (propagated with 
different coefficients of absorption) are linearly polarised along the principal 
planes of absorption. Similarly, the waves propagated along directions 
appreciably inclined to the optic axes will be practically plane polarised along 
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the principal planes of refraction, since the absolute values of the dichroism 
(k,—k,) is usually such that it is very small compared with the birefringence 
(5,;— 8) along such directions. Also, where the principal planes of absorp- 
tion and refraction coincide (as for example along the axial plane in ortho- 


rhombic crystals) the waves will be rigorously linearly polarised along the 
common principal planes. 


The more general case, where the solution is not so apparent, is dis- 
cussed analytically in Section 7, but the main results will first be proved 
more briefly and elegantly by the use of the Poincaré sphere. For this 
purpose, the form of the arbitrary elliptic vibration P—as distinct from its 
intensity and absolute phase—must first be specified by means of certain 
parameters. The principal planes of absorption and refraction form the 
two natural co-ordinate systems to which the vibration P may be referred. 
The ratio tan ¢ of the amplitudes of the components of the vibration P along 
OY, and OX, does not by itself completely specify the form of the elliptic 
vibration (since the phase difference 6, between these components has also 
to be given). Similarly the ratio tan % of the amplitudes of the components 
of P along OY; and OX, does not by itself completely specify the form of 
P (since the phase difference 0, between these components has also to be 
given). But ¢ and ¢% together form two convenient symmetrical parameters 
completely specifying the form oi the elliptic vibration P (¢, %)—provided 
we separately give the sense of description of the ellipse. It may also be 
noted that apart from an intensity factor, cos 2¢ and cos 2% are the two 
values of the second Stokes parameter M of the vibration P (4, +) when it is 
referred successively to co-ordinate systems along the principal planes of 
refraction and absorption respectively. 


As indicated in the figure, it turns out that there are two particular 
elliptic vibrations described in the same sense, Pg (¢q, a) and Pp» ($5, ¥p), that 
can be propagated without change of form under the superposed effects of 
birefringence and dichroism. The form of the vibration Py can be obtained 
from that of Pg merely by rotating the latter by 90° in its own plane—which 
means that 4p and %p are complementary to ¢q and wg respectively. 


4. Use OF THE POINCARE SPHERE FOR SUPERPOSITION 


(a) The General Method 


The Poincaré sphere,“ which has proved very useful for the analysis 
of the propagation of polarised light in ‘ransparent media, turns out also to 
be Of great use in our present discussion on absorbing crystals. 
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As is well known, a one-to-one correspondence can be set up between 
all the points on the surface of a sphere (the Poincaré sphere) and all the 
possible forms of elliptic vibrations that can be conceived (circular and linear 
vibrations being regarded as particular cases of elliptic vibrations). In 
particular the arbitrary elliptic vibration P referred to in the previous section 
is represented by a corresponding point P on the Poincaré sphere, while a 
linear vibration along OX; will be represented by some other point X;. The 
infinitesimal operation of anisotropic absorption described in the last section 
will obviously alter the form of the elliptic vibration P in such a manner 
that one may say it gets more polarised in the direction of OX;,, since this is 
the less absorbed component. This infinitesimal alteration in the state of 
the elliptic vibration P corresponds (on the Poincaré sphere) to an infinite- 
simal movement of the point P directly towards Xx, i.e., along the direction 
of the shortest arc joining them. 


Similarly the infinitesimal alteration in the state of the (initial) elliptic 
vibration P due to the operation of birefringence alone, corresponds to an 
infinitesimal movement ds, of the representative point P. If the elliptic 
vibration is to be propagated without change of form, this movement ds, 
should be equal and opposite to the displacement of P due to dichroism 
alone; and the problem of finding the states of polarisation that can be 
propagated without change of form is therefore reduced to the simple geo- 
metrical problem of finding the points P on the Poincaré sphere which 
satisfy the above requirement. 


(b) The Operations of Dichroism and Birefringence 


Referring to Fig. 2, let X,_ and Y, give the orientations of the principal 
planes of absorption. (The arbitrary elliptic vibration P has not been indi- 
cated on the sphere, but for the sake of concreteness, may be temporarily 
identified with the particular state Pg in the figure.) If the elliptic vibration 
P is resolved into two orthogonal linear vibrations in the states X; and Yj, 
then the amplitudes F;, and Gx of these components will be proportional 
to cos % and sin %, where 2% is the angular distance of the point P from X, 
on the Poincaré sphere. (For a proof of this statement, see reference 7.) 
Hence G;/F;, = tan %. If the amplitudes of these components are reduced 
by the multiplying factors e* and e-**, the angular distance of P from X; 
will change from 2% to 2x’ where: 

— kz 
tan y’ = BEE = tan pethrk? (1) 





Since the phases of the X; and Y,; components of P are to be left unaltered 
in this operation, the movement of P will be entirely on the meridional arc 
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Fic. 2 
(X,, Y,) and (X,, Y,)—Principal planes of absorption and refraction, respectively. 
P, and P, —States of polarisation propagated unchanged. 
X,X_ = 2X; XyX’ = 2X; XX’ = 2X,; PX’ = 28. 


Y;,PX;. This follows from the fact that the X;, and Y;, components of all 
elliptic vibrations on this arc have the same phase difference 6), (where 9, is 
the angle indicated in Fig. 2); for, if an additional phase difference — 4, 
be introduced between these components, any such elliptic vibration will 
be reduced to a linear vibration on the equatorial arc X;Y,;Y; (by a well 
known property of the Poincaré sphere). 


The infinitesimal operation of linear dichroism (corresponding to a 
passage dz) will, apart from reducing the intensity, cause the initial state of 
polarisation P to move along the arc PX; towards the state X; (the less 
absorbed component), through an arc ds, = — 2d. From (1) we have 


tan % + d(tan #) = [1 — (k, — k,) dz] tand 


or 


2 sec? pdb = — 2(k, — k,) dz.tan % 
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which leads to the simple relation 


dsy. = (k sin 2h) dz (2) 


where 
k denotes (k, — k,) 


Let X, and Y; represent the orientations of the principal planes of 
refraction, the former corresponding to the slower wave in the absence of 
absorption; let 5 denote the phase difference (5, — 6.) introduced per unit 
distance (in the absence of dichroism). The infinitesimal operation of bire- 
fring:nce (corresponding to a passage dz) consists in rotating the sphere 
clockwise about the equatorial diameter X, Y; through the infinitesimal 
angle 5dz. This operation will cause the initial state of polarization P to 
move along the arc of a small circle with X; as pole through an arc ds,, where 


ds, = (6 sin 2¢) dz (3) 


(c) The States of Polarisation Propagated Unchanged 


In order that the simultaneous superposition of linear dichroism and 
birefringence should cause no change in the state of P, the movements dsj, 
and ds, must be equal in magnitude and opposite in direction. Since arc 
ds;, is along PX, while arc ds, is perpendicular to PX; we must have firstly, 
X, PX; = w/2, or 

cos 2X = cos 2¢ cos 2% (4) 
and secondly 
5 sin 2¢ = k sin 2 (5) 


together with the condition that P will be a right- or left-elliptic vibration 
according as X (the angle between the X, and X; axes) is positive (0 to 7/2) 
or negative (0 to —7/2). 


In general there are two positions P which simultaneously satisfy these 
conditions—the relations (4) and (5) being unchanged when we alter 2¢ and 
2 to 7 — 2¢ and 7 — 2% respectively. Thus in the figure the state Pp whose 
distances from the points Y; and Y; are 2% and 2¢, is also propagated un- 
changed. The states Pg and Pp have the same latitudes, their longitudes 
differing by 7. 


Hence we arrive at the result, also obtained from the electromagnetic 
theory, that the states of polarisation propagated unchanged along any 
general direction are two similarly rotating elliptic vibrations which have 
their major axes crossed and their ellipticities equal, 
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In order to construct these two elliptic vibrations (see Fig. 1) we must 
first determine the orientations of one of the principal diameters—for 
example the orientation of OX’ which is the major axis of Pg and also the 
minor axis of Pp. Let the inclination of OX’ be X, (anti-clockwise) with res- 
pect to OX, and X, (clockwise) with respect to OX;,. The direction of OX’ 
may be determined by the relation 


sin 4X, ___ 8? 
sin 4X, k? (6) 


Next the ratio tan @ of the minor to the major axis of the elliptic vibrations 
may be obtained from 


sin? 20 = tan 2X, tan 2X, (7) 


—a relation which gives the ellipticity (tan 9 = cosec*20 + 1/cosec? 20— 1) 
in terms of the orientation of the axes of the elliptic vibrations. It may be 
noted that equations (6) and (7) remain unaltered when we change X, and X, 
to (7/2 + X,) and (7/2 + X,) respectively. 


The relations stated in the previous paragraph will now be proved with 
the aid of the Poincaré sphere (Fig. 2). The direction OX’ is given by the 
point X’ on the equator having the same longitude as P. Then since the 
triangles X,PX’ and X;.PX’ are both right-angled, we have 


cos 2¢ = cos 2X, cos 20 
cos 2s = cos 2X, cos 20 J (8) 


Multiplying these equations and comparing with (4), 
cos (2X, + 2X5) = cos 2X, cos 2X, cos?20 
which on simplification gives the relation (7). 
To prove (6) we consider the right-angled triangles X,PX, and X;,PX’: 


- tan 2 _ tan 2X, 


cos X_ = tan 2X tan 2% 
Hence 

tan? 2.4 = tan 2X, tan 2X 

tan? 2¢ = tan 2X, tan 2X 
or, 


tan? 2% _ tan 2X, 9 
tan? 2¢ tan 2X, (9) 
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From relations (8) we also have 


cos? 2% _ cos” 2X, 
cos? 26 ~ cos? 2X, (10) 
Multiplying (9) and (10) and comparing with (5) we get the required rela- 
tion (6). 


(d) Comparison with the Electromagnetic Theory 


The relation (7) giving the ellipticity in terms of the orientation of the 
axes of the elliptic vibration, is identical with that deduced from the electro- 
magnetic theory (Pockels, Joc. cit., p. 399, eq. 54); while the relation (6) 
giving the orientation of the axes of the elliptic vibrations has to be com- 
pared with the following similar relation (Pockels, loc. cit., p. 397, eq. 53): 


sin 4X, __ p? ' 

sin 4X,  ? (6’) 
where 

p =4(d, — a,) and o = $(b, — dy) 


The tensor components 4;, da, etc., may be easily shown to have the 
following geometrical meanings. The major and minor semi-axes of the 
elliptical section of the index ellipsoid made by the plane of the paper have 
lengths 1/4/a, and 1/+/d, respectively, while the major and minor semi-axes 
of the elliptical section of the absorption ellipsoid have lengths 1/./b, and 
1/,/b, respectively. Relations (6) and (6’) will be identical if 


8°/k® = p®/o® (11) 


As pointed out in Section 2, the waves propagated along directions 
appreciably inclined to the optic axes may be considered as linearly polarised ; 
and for such directions of propagation, if 1/./b be the length of any radius 
of the absorption ellipsoid, then b = 2«v3/c, where « is the extinction co- 
efficient and v the velocity for that vibration-direction. In consonance 
with this it would be natural to use the following relation for the hypothetical 
extinction coefficients «, and «x, in the absence of birefringence: 


3 3 
21m wiles b, : 2k Vm Si bs (12) 


c c 
where vm is a mean velocity. We then have 


2a 2mc o 
k= o (Ky — Ky) = % Tan? 
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Relation (11) will be obviously satisfied if we analogously set 


_ 2ac pP _ 2m¢ a@—a 
">. a. a (13) 


Since the velocoties v, and v, in the absence of absorption are equal to \/a, 
and 1/d. respectively, relation (13) will be exactly satisfied if we define the 
mean velocity vm by 


Vm> = $ (Yy + Vq)°VyVe (14) 


5. THE ABSORPTION COEFFICIENTS AND REFRACTIVE INDICES OF THE WAVES 


It is well known that in the case of a transparent crystal, it is simpler 
to specify the velocities of the waves as functions of the vibration directions 
than as functions of the directions of propagation: the former leads to the 
simple index-ellipsoid representation, the latter to the comparatively more 
complex wave surface of two sheets. We shall show that for absorbing 
crystals too, if we choose to express the velocities and absorption coefficients 
as functions of the states of polarisation (4, ¥) of the waves, the resulting 
expressions (as deduced both by the method of superposition and by the 
electromagnetic theory) may be put in a very simple form. 


When an elliptic vibration of unit intensity in any state of polarisation 
travels a distance dz, the diminution in its intensity may be calculated directly 
from the reduction of intensity involved in the infinitesimal operation of 
dichroism corresponding to the passage dz—since the operation of birefring- 
ence produces no reduction in intensity. If, in addition, the elliptic vibration 
be in astate of polarisation Pg that can be propagated without change 
of form, this reduction in intensity may be equated to 2k, dz where kg is 
the coefficient of absorption for that wave. The amplitudes of the X; and 
Y;, components of the elliptic vibration Py will be cos %q and sin yg respec- 
tively; hence the reduction in intensity of these components will obviously 
be 2k, cos? %q dz and 2k, sin* aq dz respectively. Therefore, 

ka = k, cos*hq + kz sing 
Similarly | 
kp = k, sin® tha + ke cos? tq (15) 
So that 
(ka — ky) = (ky — ke) cos 2 


Here 2); being the arc PX, on the Poincare sphere may be evaluated by 
the relations (4) and (5) of Section 4c, which determine the states of polarisa- 
tion of the waves, 
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Expressions for the refractive indices of the waves (in terms of the state 
of polarisation of one of them) are equally simple, being given by: 


Ng = N, COS* dq + Ne sin? da 
Ny = N, sin? dg + Ny cos? da (16) 
(na — No) = (m — Nz) Cos 2¢q 


The proofs of these relations do not have the same simplicity as those 
of (15), and will be given only at the end of Section 7, since recourse must 
be taken to the analytically derived equations obtained there. 


(6) THE PROPAGATION OF LIGHT ALONG THE AXES OF CIRCULAR 
POLARISATION 


(a) The Singular Axes 


The electromagnetic theory predicts that close to an optic axis and on 
either side of it, there exist two directions along each of which only one state 
of polarisation (and not two) can be propagated unchanged: only a right- 
circularly polarised wave can be propagated along one of these axes, and 
only a left-circularly polarised wave along the other. These directions have 
been termed the Windungsachsen; Voigt has also referred’ to them as 
singular axes and we shall follow this simpler nomenclature. At these two 
axes of circular polarisation the inclinations of the principal planes of absorp- 
tion with respect to the corresponding principal planes of refraction are +45° 
and —45° respectively. Further, along these two directions the pure bi- 
refringence term 5 is equal to the dichroic term k. 


The remarkable property of these axes follows very simply from the 
standpoint of the method of superposition by the use of the Poincaré sphere. 
Let us suppose for example that the principal plane of absorption OX;, makes 
an angle of —45° with respect to the corresponding principal plane of refraction 
OX,. In Fig. 3, the diameter X,Y, will then be at right angles to X,Y; as 
shown. If we consider a state of polarisation initially coincident with the 
pole Cj, it can be seen that its movement ds, (due to an infinitesimal clock- 
wise rotation 5dz, about X,Y,) will be oppositely directed to the movement 
ds), towards the less absorbed component X;; and the movements will be 
equal in magnitude if =k. Thus a left-circular vibration can be propa- 
gated unchanged along such a direction. Further, there can be no other 
state which can also be propagated unchanged, since the 2 elliptic vibrations 
propagated unchanged along any direction must have the same sense of 
description and the same ellipticity. Similarly, where OX; makes an angle 
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Ce 











Xp 





Cr 
Fic. 3 
C,—Right circular vibration incident in the direction of a singular axis (where only the 
left-circuiar vibration C, can be propagaied unchanged). 
Q—The state of vibration at a depth z (specified by the arc s). 


of + 45° with respect to OX;,, and where in addition 5 = k, only a right- 
circular vibration can be propagated unchanged. The refractive index of 
the circularly polarised wave that can be propagated unchanged along a 
singular axis is 4(m,-+ mg) and its absorption coefficient 4(k, + k,) as 
may be seen by setting ¢ = 4 = 7/4 in relation (18) and (19). 


Before proceeding to discuss in more detail the propagation of light 
along the singular axes, we consider it relevant to point out that the func- 
tions with which we are concerned show no discontinuity at the singular 
axes. Thus both the elliptic vibrations propagated without change of form 
along any general direction, gradually degenerate into two (identical) circular 
vibrations as we approach a singular axis from any side whatsoever. [This 
can be seen by making X > 45 and 5k in relations (4) and (5) of Sec- 
tion 4c.] The refractive indices and absorption coefficients of these two 
waves, being determined by their states of polarisation (by relations 15 and 
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16 of Section 5) tend towards the common values 4 (nm, + m2) and 4 (k, + kg) 
respectively, as we approach a singular axis. (See also reference 8.) 


(b) Effects with Incident Circularly Polarised Light 


In this section we shall inquire as to what will happen when, for example 
right circularly polarised light C, is incident in the direction of a singular 
axis where only a left circular vibration C; can be propagated unchanged 
(Fig. 3). Our results in this connection are at variance with those expected 
by Voigt. It was supposed by Voigt* * that if a plate cut normal to an optic 
axis is viewed in convergent circularly polarised light, then along the singular 
axis where the incident vibration can be propagated unchanged, more light 
would get through than in the neighbourhood of the other singular axis 
where only the oppositely directed circular vibration can be propagated 
unchanged; and that the latter direction should in consequence appear 
darker than the former.* On performing an actual experiment, he observed a 
dark and a bright spot in the field of view, one on either side of the optic 
axis and this was considered by him as confirming his view. According to 
our analysis, however, it is the singular axis where the incident vibration 
can be propagated unchanged that should appear darker than the other 
singular axis (where only the oppositely directed circular vibration can be 
propagated unchanged). 


We shall apply directly the method of superposition, according to which, 
given the state of vibration at a particular plane in the medium, the state of 
the vibration at a further distance dz is obtained by superposing the effects 
of pure birefringence and pure dichroism corresponding to that passage. 
The state of vibration should then get progressively altered as we proceed 
into the medium. 


Referring to Fig. 3, if the state of polarisation be initially coincident 
with the pole C,, its movement ds, (due to a clockwise rotation 5dz about 
X,Y;) is in the same direction as its movement ds; towards the less absorbed 
component X;; and the sum of these movements will give the alteration 
in the state of vibration corresponding to a passage dz. Continuing this 
procedure, it can be seen that as we proceed into the medium, the state of 
polarisation progressively moves along the arc C,;X;C;. At a particular 
depth the vibration would be linearly polarised along the principal plane 


* In a later paper ®, Voigt has suggested that if we could get a plate exactly norma! to a singular 
axis, and have circularly polarised light of the proper sense incident precisely along this normal, 
the light would be totally refiected—the reflection being partial in practical cases. This idea receives 
no support from the results of the present investigation. 
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of absorption OX,. After this stage the movement ds; due to dichroism 
opposes the movement ds, due to birefringence; but the latter being greater 
in magnitude, the state of polarisation continues to alter as we proceed 


further into the medium, tending towards the state C; that can be propagated 
without change of form. 


Thus if we consider the state of polarisation at successive depths within 
the medium, we see that the incident right-circular vibration will first get 
modified to an elliptic vibration (with major axis always at 45° to the principal 
planes of refraction), which in turn gets reduced to a linear vibration; as 
we proceed further the linear vibration opens out into a /eft-handed elliptic 
vibration, which gradually tends towards the state of a left-circular vibration 
that can be propagated unchanged. Nevertheless, as may be seen physically, 
this last state is never attained at any finite depth; for as the state of polarisa- 
tion comes close to that of a left-circular vibration, the modification of the 


state corresponding to an additional passage dz becomes correspondingly 
reduced. 


It is easy to deduce an explicit expression for the state of polarisation P 
that should be expected (according to the above line of argument) at any 
depth z inside the medium. The state P may be specified by giving the length 
s of the arc C;P. Then the state s + ds at the depth z + dz will be given by 


ds = (8+ k cos s) dz (17) 


according to relations (2) and (3) of Section 45. Since 5 = k, we have on 
integration, 
tan 4s = kz (18) 


This relation shows that the transformation from a right-circular vibration 
(s = 0), to a linear vibration at 45° to the principal planes of refraction 
(s = 7/2), occurs within a smaller depth than if the crystal had been trans- 
parent; whereas the corresponding alteration from the linear vibration to 
a left circular vibration (which for a transparent crystal would have occurred 


at a finite depth) requires here an infinite passage, due to the ‘ retarding’ 
effect of the dichroism. 


We shall next calculate the intensity I, of the vibration P at a distance z 
inside the medium. The diminution of intensity — dIz corresponding to 
an additional passage dz is given by 

— dI-/I- = 2k-dz (19) 


where, it must be noted k, is not a constant but a function of the state of 
polarisation and hence also of the depth z. We will have for kz an expression 
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analogous to (15), Section 5: 
kz = k, cos? # + ky sin? % 
= $(k, + ke) — tk cos 2% (20) 


Thus kz is always less than the coefficient of absorption 4 (k, +k.) of 
the left-circularly polarised wave that can be propagated without change of 
form along the same direction. (This is more directly seen by the fact that 
the state of polarisation is always nearer to the less absorbed component X, 
than is a left-circular vibration.) Hence when the sense of description of a 
circular vibration incident in the direction of a singular axis is opposed to that 
which can be propagated unchanged along that direction, the emergent intensity 
should in fact be greater than when the sense of description of the incident 
vibration is reversed. An expression for the ratio of the emergent intensities 
in the two cases will now be deduced. Since, from our point of view, the 
incident disturbance can propagate into the medium in both cases (though 
in one case with a progressive change in the state of polarisation) we have 
no particular reason to assume that the reflection losses would be different 
in the two cases. 


Substituting the value of kz given by (20) in (19) we get 
— diz/I, = (k, + k.) dz — k sin s dz 


Expressing sin sin terms of z by using relation (18), and integrating, we have 
—if 1, be the emergent intensity and I, the intensity entering the medium, 


log (Io/I1) = (ky + Ke) z — log (1 + k?z*) (21) 
On the other hand if I, be the emergent intensity when the incident circular 
vibration is of the sense which can be propagated unchanged, 

log (Io/I2) = (ky + ke) z (22) 


From (21) and (22) we have the following simple relation for the ratio of the 
intensities emerging in the two cases: 


I, /I, = 1 che k?z? 
a ratio which is always greater than unity. 


Though our results regarding the properties of the singular axes are at 
variance with those expected by Voigt, it must not therefore be concluded, 
that the method of superposition leads to results differing from the electro- 
magnetic theory—since it is possible to regard the former merely as a mathe- 
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matical artifice: for (as has been pointed out in Section 8), by reversing the 
entire scheme of arguments, the propagation of two elliptically polarised 
waves with different absorption coefficients and velocities (as given by the 
electromagnetic theory) may for mathematical convenience be rigorously 
treated by the method of superposition; the convenience of this artifice 
being particularly manifest when we wish to find the limiting effects as the 
two elliptically polarised waves gradually tend towards the state of two 
(identical) circularly polarised waves with equal velocities and absorption 
coefficients. 


7. ANALYTICAL DISCUSSION OF SUPERPOSITION 


Referring to Fig. 1, let us suppose that we are given the equation of the 
elliptic vibration described at any particular plane z in the medium. [If the 
initial state of polarisation is to be propagated without change of form, then 
the equation of the vibration at the plane z + dz can be obtained not only 
by the method of superposition but also from the usual equation for the 
propagation of a damped wave. By equating these two expressions we can 
determine not only the states of polarisation that can be propagated without 
change of form but also their velocities and extinction coefficients. 


Let OX and OY be two arbitrary rectangular axes taken in the plane 
of the figure, the inclinations of OX; and OX, with respect to the positive 
x-axis being a, and o, respectively. Let the components of the arbitrary 
elliptic vibration P along the axes OX, OY, have the following equations 
(using complex notation and indicating the complex quantities by bars): 


x = Fellt-o) — feirt 
y= Gei(“t—62) — get (23) 
so that 
% aN S gia 


where 2/ f is the ratio of the complex amplitudes, G/F the ratio of the real 
amplitudes and @ the difference of phase (0, — 0.) between the x and y 
components. (It may be noted that G/F and @ will have—on the Poincaré 
sphere—geometrical interpretations essentially similar to those that have 
been described in Section 4 c for G;,/F, and 9%.) 

Let the complex amplitudes become /’, &’, after the infinitesimal ope- 
ration of dichroism alone, and /”, 2”, after both the infinitesimal operations 
of dichroism and birefringence corresponding to a passage dz. But if the 
initial state of polarisation (/, Z) is to be propagated without change of form, 








o 
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and with a specific extinction coefficient « and refractive index n, then its 
state (f", 2”) after propagating a distance dz should also be given by: 


frei = fe (wt — ~ idz) 
f =f (1-1 2 ad) 


and 
1 ns » 
‘= 1-—i;aA dz ) i 
# =e (1-i< 
where 7 is the complex refractive index (n — ix) of the elliptically polarised 
wave. We shall now, by the method of superposition, proceed to determine 


expressions for the final state (7, 2”) in terms of the initial state ( f, 2), and 
then substitute these expressions in (24). 


The Operation of Dichroism.—The elliptic vibration P, given by (23), 
is first referred to the axes OX;,, OY; along the principal planes of absorp- 
tion; the complex amplitudes (f;, 2,) of the components along these direc- 
tions can be obtained from the amplitudes (7, Z) by the usual transforma- 
tion scheme for the rotation of co-ordinate axes through an angle a,. These 
amplitudes 7; , 2, are then multiplied by (1 — 2m/Ay. «,dz) and (1 — 2m/Ag- 
k,dz) respectively, to give the amplitudes f;,’, 2,’ of the XK, Ye 
components after the operation of dichroism. Finally the elliptic vibration 
Jx', Z;’ thus obtained is referred back to the axes OX, OY; the complex 
amplitudes 7’, 2’ of the components along these directions are obtained from 
the amplitudes 7;,’, 2,’ by the usual transformation scheme for the rotation 
of the co-ordinate axes through an angle — a,. 


If we go through the straightforward but lengthy algebraic substitutions 
involved in the procedure described in the last paragraph, it can be shown 
that the amplitudes 7’, 2’ of the x, y components after the operation of 
dichroism are related to the corresponding initial amplitudes 7, 2 by relations 
which may be put in a form analogous to (24): 


at (1% (ou + fone 


or 


(24) 


. 7 (25) 
os = G5 _ — uf 
& =g8 [1 p™ (ze oo Fa Kis) dz | 
where, if a, be the orientation of the positive OX, axis, 
Ky, = K, COS? a, + xg sin? a, t 
(26) 


Keg = K, Sin? a, + kK, COS? a, 


Ki9 — 5 (ky —_ Ka) sin 2a, 
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If an ellipse be drawn with its principal semi-axes of lengths 1/,/x, 
and 1/+/«. lying along OX; and OY,, then 1/4/«,; and 1/+/«g, are the lengths 
of the radii vectors intercepted by the OX and OY directions, the equation 
of the ellipse being 

Ky X? + Kaa Y? + 2ky2 XY = | 

The Operation of Birefringence-—The mathematical procedure involved 
in the infinitesimal operation of birefringence is essentially the same as in 
the operation of dichroism. The elliptic vibration /’, Z’ obtained after 
the operation of dichroism is first referred to the axes OX,, OY, along the 
principal planes of refraction; the complex amplitudes /,’ Z,’ of the compo- 
nents along these directions are then multiplied by exp (— i 2m/A).n,dz) 
and exp (— i27/Ay.n,dz) respectively—where n, and nz are the refractive 
indices in the absence of absorption. On referring the final vibration 
(fr", Zr’) back to the axes OX, OY we will have the complex amplitudes 
(f", 2") of the x, y components, related to the corresponding amplitudes 
f',2' (before the birefringence operation) by equations essentially similar 
to (25), though put in the form of (24): 


f'=f' [1 = i< (mu + & . my.) dz | 
2” =2' [1 -i< (m+ & my) dz | 


Here 74, M22 and M4. are to be regarded as defined by the relations (ana- 
logous to (26): 


(27) 


Ny, = Ny, COS? ag + Nz Sin? ay 
Nog = Ny Sin® ag + Nz COS? ay (28) 
Ny, = 4 (m, — ng) sin 2a, 

where a, gives the angle made by the positive OX, axis with the x-axis. 


Since we shall omit terms involving dz*, the value of (2’/f’) to be substi- 
tuted in (27) need not include even the terms of the first order in dz, i.e., we 
may write (2/7 ) for (2’/f ’) in (27). We then obtain as the equation connecting 
the final state of polarisation f”, 2” (after both the infinitesimal operation of 
birefringence and dichroism corresponding to a passage dz), with the initial 
state f, Z: 


| a 


I 


F{1 13 [na — ten + § + ta teas | de 
A (29) 


¥ g{1-i< [a = een + 7 * tha — teas | dz } 
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The Final Equations for Wave Propagation——We now introduce the 
values of 7”,2” given by (29), into equation (7). Conciseness will obvi- 
ously be attained if we first introduce the complex quantities: 

Ay = My — Iky 

Nog = Nag — ike 

Tyg = Nyg — Ikyg 
where the mx and «px have already been defined in relations (26) and (28). 
We then obtain as our final equations: 


7 (30) 


PT) 
| 
= 
ie 
I 
' 
=! 
ey 


There will be two pairs of values (nq, Za/fa) and (ip, Zp/f,p) which simulta- 
neously satisfy (30); and—since Z/f and fi are both in general complex— 
this means that there should be two elliptically polarized waves that can 
be propagated, each with a specific velocity and coefficient of extinction. 
Eliminating 7 between the two equations of (30) by subtracting, we get the 
following quadratic in g/f determining the states of polarisation propagated 
without change of form: 

ft g = Nyy as Ng 
of G) 

The two roots of this equation are obviously connected by the relation 
(Zalfa) = — (fv/%d), from which it follows that the two elliptically polarised 
vibrations have their major axes at right angles and their ellipticities equal, 
but are described in the same sense (see, e.g., McLaurin’). 


Eliminating 2/f between the two equations (30) by multiplying the two, 
we get the following quadratic in # determining the complex refractive 
indices of the waves: 


(fi — fy) (A — figg) = 1,2” (32) 


From this we get the expressions for the sum and difference of the complex 
refractive indices: 


(fig — fp)? = (finn — fgg)? + 4iy2" (33) 
and (lig + Np) = Ty + Nog (34) 
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Velocities and Absorption Coefficients of the Waves.—The expressions for 
the velocities and absorption coefficients given in Section 4 may be derived 
from the equations (30) which give the complex refractive 7 in terms of the 
corresponding state of polarisation Z/f, If we choose axes of co-ordinates 
along OX;, OY, then 7, = 0, ag =0 and a, = 2X. The first of the two 
relations in (30) gives: 


(nt — ie) = (1%, — iy) + SF e¥@r - ( — ine) 
T 


Equating real parts 
n =n, + (G;/F,) ky, sin Or 
= n, + (G,/F,) - 4 («, — «2) sin 2X sin 6, 
Referring to Fig. 2, since 0, = PX, X; we have 
sin 2X sin 0, = sin 2% 
Hence on using eq. (5), Section 4c, we have 
n =n, — (G;/F,) - $(m% — ng) sin 26 


Since (G,/F,) = tan ¢, we get as our final expression for the refractive 
index, 


= n, cos? ¢ + ng sin? d 
Similarly we will have for the extinction coefficient, 
kK = kK, Cos? f + Kz sin? o 
8. COMPARISON WITH THE ELECTROMAGNETIC THEORY 


Let the sections of the index and absorption ellipsoids made by the 
xy plane be given by the respective equations: 


Qy,X* + Agey® + Zdyxy = 1 35 
basx* + Dysy? + Dua xy = I " 
Let us introduce the quantities 

Chk = Gnk + ibn (36) 


Then the equations (30) giving the states of polarisation and the complex 
refractive indices of the two waves propagated along any direction have to 
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be compared with the following similar relations obtained from the electro- 
magnetic theory: 


(v? — ew) = 3 Cie 


(37) 
(v? — Cea) =5 * Cra 


where the complex velocity v = v (1 + ix’) 


We may first remark that it would indeed be possible to define the quanti- 
ties 71,1, figg and fiyg (in terms of ¢,, Cog and Cy) in such a manner that the 
results obtained by the method’ of superposition would be identical with 
the results of the electromagnetic theory. But in order to retain the physical 
content of the method of superposition it is necessary to regard the velocities 
v, and v, in the absence of absorption as being equal to 1/a, and +/dg respec- 
tively (where 1/4/a, and 1/+/a, are the lengths of the principal radii of the 
elliptical section of the index ellipsoid). And once this is done, at least some 
of the results obtained by the method of superposition have necessarily 
to be regarded as approximations. We shall however show that for direc- 
tions near an optic axis where the birefringence is necessarily very small and 
where alone the ellipticity of the waves play an important role, the error 
involved is negligible. To this end we shall start by assuming the relations 
(12) and (14) which give a connection between the extinction coefficients 
k, and x, on the one hand, and the lengths 1/+/b, and 1/+/b, of the principal 
radii of the elliptical section of the absorption ellipsoid on the other. The 
quantities an, occurring in equation (35) can obviously be expressed in terms 
of a, and a, thus: 


ay = ay, cos? Qo + ag sin? Ag 
Ag = A, Sin” a, + dy COS? ay (38) 
ayo = 4 (a, reas as) sin 2a 


Similar relations analogous to (26) hold for the bp x. 


On examining the equations [(30) and 37)] obtained by the method of 
superposition and by the electromagnetic theory, we notice that the two are 
entirely similar in form, the only difference being the occurrence of the 
quantities cpz instead of ip~, and v* instead of 7. Hence it follows that given 
any equation obtained by the method of superposition, a corresponding 
exact equation obtainable from the electromagnetic theory can be written 
down, merely by changing the symbols occurring in the equation according 
to the following scheme; 
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n ny Ng K Ky Ko 
3 
2 2xv 
° a ay saaitie ~ —b, —b, 


Thus to obtain the states of polarisation (¢, 4) that can be propagated with- 
out change of form, we have to replace the ratio 5/k in eqn. (5) by the ratio 
plo where p = 4 (a, — a,) and o = 4(b, — b,). This replacement will how- 
ever leave the equation unaltered, as we have already shown in Section 4 (d) 
that 5/k = p/o. It is also possible to show that the sense of description of 
the two vibrations as obtained by the method of superposition is the same 
as that obtained by the electromagnetic theory. 


Expressions for the velocities and absorption coefficients according to 


the electromagnetic theory may be similarly written down from the relations 
(15) and (16): 


v? = a, cos? ¢ + a, sin® (39) 
3 
240" = by cos? yh + by sin® y (40) 


We hence obtain 
(Yq? — Vp”) = (a, — ag) COS 2¢q 


Or, if Ym be the mean velocity introduced in Section 4 (d) 


Vp? — Yq? 
‘—— “= (n, — Ng) cos 2¢ 


On comparing this with the last equation in (16) we see that the approxima- 
tion involved in using the method of superposition is to regard the expression 
on the left-hand side of the above equation as being practically equal to the 
difference in refractive indices (mq — np). This is justifiable along directions 
where the birefringence is low, and in fact this same approximation is also 
made when the propagation in transparent optically active crystals is re- 
garded from the standpoint of superposition (Pockels, /oc. cit., p. 312). 


From Egn. (40) we get 


2kaa? — 2xpvy®? = (b, — bg) cos 2hg 
Or, using (12) 


Kqlq* —_ Kpup® 


0.3 = (ky — Kg) cos 2a 
™ 


On comparing this with the last equation in.(15), we see that the approxima- 
tion we have to make is to regard the expression on the left as being prac- 
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tically equal to the difference in extinction coefficients («q — <p)—an approxi- 
mation that is again justifiable where the birefringence is low. 


It may be similarly shown that the mean refractive index and the mean 
extinction coefficient are also only negligibly in error. We shall not how- 
ever give the proof here since the absolute values of the velocities and extinc- 
tion coefficients are not as important as their differences. 


The author wishes to express his thanks to Professor Sir C. V. 
Raman for his kind interest in this work. 


9. SUMMARY 


The optical behaviour of pleochroic biaxial crystals in the vicinity of 
an optic axis can be elegantly interpreted as due to the effects of linear 
birefringence and linear dichroism superposed continuously along the depth 
of the material. This idea—followed up geometrically with the Poincaré 
sphere, and analytically by direct algebraic methods—explains the 
elliptical polarisation of the two waves propagated in any general direction, 
as due to the non-coincidence of the principal planes for the usual operation 
of birefringence with those for the operation of dichroism. In particular, 
it also explains the existence, on either side of an optic axis, of two singular 
axes—along any one of which only one circular vibration with a definite 
sense of description can be propagated unchanged: for these directions 
the principal planes of absorption and refraction make angles of 45° with 
each other, the linear birefringence and dichroism being also equal in 
magnitude. 


A discussion is also given of the effects to be expected when a circular 
vibration incident in the direction of a singular axis has a sense of descrip- 
tion opposite to that which can be propagated unchanged; and it is shown 
that—contrary to what Voigt expected—the emergent intensity will not be 
either zero or negligible, but will in fact be greater than when the sense of , 
description of the incident vibration is reversed. 
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